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1 Page  6,  Equation  2 should  read 

N 

v(C»n)=B  +n+y'  (B  cosh  nkn  + C sinh  nkn)  cos  nkC 
o ",  n n 

n=l 

2.  Page  19,  lines  2 and  3 should  read  "and  Fjj  are  the  point  values  of 

F from  Equation  11  at  Qs  , Qf  > H points,  respectively, 

is  the  point  value  of  p from" 

3.  Page  28,  the  second  line  following  Equation  k9  should  read  "a  is  a 
rotation  angle  used,  primarily,  to  fine  tune  the  alignment  of" 

Page  i*5,  line  22  should  read  "raately  0900Z,  8 September  (between  the 
time  of  Plates  70  and  7l)-  The" 
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The  investigation  reported  herein  was  conducted  under  Department 
of  the  Army  Project  No.  i*Al6ll01A91D,  In-House  Laboratory  Independent 
Research  (ILIR)  Program.  The  program  is  sponsored  hy  the  Assistant 
Secretary  of  the  Army  (R&D).  The  research  was  performed  at  the  U.  S. 

Army  Engineer  Waterways  Experiment  Station  (WES),  Vicksburg,  Mississippi. 

The  study  was  conducted  by  Dr.  John  J.  Wanstrath  of  the  Coastal 
Branch,  under  the  general  supervision  of  Mr.  H.  B.  Simmons,  Chief  of 

the  Hydraulics  Laboratory,  Dr.  R.  W.  Whalin,  Chief  of  the  Wave  Dynamics  ' 

i 

Division,  and  Dr.  C.  L.  Vincent,  Chief  of  the  Coastal  Branch.  Dr.  D.  L.  ! 

Durham,  CAPT  F.  C.  Perry,  and  Mr.  J.  W.  McCoy  assisted  in  the  operation  ! 

and  data  collection  of  the  Murrells  Inlet  physical  model  surge  experi- 
ments. Special  assistance  in  providing  data  and  other  services  was 
given  by  the  U.  S.  Army  Engineer  District,  New  Orleans,  the  U.  S.  Army 
Engineer  Division,  Lower  Mississippi  Valley,  and  the  U.  S.  Army 
Engineer  District,  Jacksonville.  Prof.  R.  0.  Reid  of  the  Department 
of  Oceanography,  Texas  A&M  University,  College  Station,  Texas,  con- 
tributed to  the  initial  formation  of  the  flooding  coast  concent  in 
open-coast  surge  computations. 

Directors  of  WES  during  the  conduct  of  this  investigation  and  the 
preparation  and  publication  of  this  report  were  COL.  G.  H.  Hilt,  CE, 

and  COL  John  L.  Cannon,  CE.  Technical  Director  was  Mr.  F.  R.  Brown.  ■ 
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CONVERSION  FACTORS,  METRIC  (Sl)  TO  U.  S.  CUSTOMARY  AND 
U.  S.  CUSTOMARY  TO  METRIC  (Sl)  UNITS  OF  MEASUREMENT 


Units  of  measurement  used  in  this  report  can  he  converted  as  follows: 


Multi pi j 


metres 

kilometres 

metres  per  second 

square  centimetres 
per  second 

millibars 


inches 


To  Obtain 


Metric  (Sl)  to  U.  S.  Customar 


3.280839 

0.5399568 

3.280839 

0.1550 

0.011*50377 


miles  (U.  S.  nautical) 
feet  per  second 
square  inches  per  second 

po\ands  per  square  inch 


U . S . Customary  to  Metric  (Sl) 


miles  (U.  S.  nautical) 

miles  (U.  S.  statute) 

square  feet 

pounds  (force)  per 
square  inch 

feet  per  second 

miles  per  hour 
(U.  S.  statute) 

degrees  (angle) 


25.1* 

0.30l*8 

1.852 

1.60931+1* 
0.09290301* 
6891*. 757 

0.301*8 

1.60931+1* 

0. 0171+5329 


millimetres 

metres 

kilometres 

kilometres 

square  metres 

pascals 


metres  per  second 
kilometres  per  hour 


radians 
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NEARSHORE  NUMERICAL  STORM  SURGE  AND  TIDAL  SIMULATION 


PART  I : INTRODUCTION 

1.  Three-dimensional  partial  differential  equations  govern  the 
motion  of  an  infinitesimal  fluid  element.  These  equations  result  from 
basic  considerations  of  mass  conservation  and  Newton's  second  law  of 
motion.  The  assumptions  involving  the  incompressibility  and  homogeneity 
of  the  water,  negligible  vertical  accelerations  of  the  fluid  parcel, 
reasonably  uniform  horizontal  flow  over  the  fluid  depth,  and  others  re- 
sult in  the  classical,  vertically  integrated  long-wave  equations  of 
motion  and  mass  conservation.  Energy  is  supplied  at  the  free  surface, 
in  general,  through  the  action  of  the  wind  and  dissipated  at  the  sea- 
bed through  friction.  The  system  of  equations  is  time-dependent,  two- 
dimensional  in  terms  of  the  horizontal  coordinates,  and  readily  in- 
tegrable  through  nvunerical  techniques.  These  equations  are  applicable 
to  the  study  of  storm  surge  generation^  on  the  continental  shelf  and 

to  the  degree  that  the  assumptions  are  valid,  nearshore  astronomical 
tide  simulation. 

2.  In  the  past,  numerical  integration  of  the  two-  dimensional  as 

well  as  the  three-dimensional  equations  of  motion  have  been  performed 

2 3 h 

using  rectilinear  grids.  Recently,  curvilinear  coordinates  ’ ’ have 
been  employed  in  two-dimensional  models,  with  particular  applications 
to  free  and  forced  long  wave  simulations  for  large  (hundreds  of  miles*) 
open-coast  stretches  of  the  continental  shelf. 

3.  Open-coast  curvilinear  models  are  considered  superior  to 
rectilinear  models  because  rectilinear  models  must  represent  the  coast 
boundary  as  a series  of  straight-line  segments  connected  at  right  angles. 
Spurious  oscillations  are  injected  into  the  calculation  by  the  boundary. 
Furthermore,  the  stair-step  boundary  must  retain  more  water  over  that 

* A table  of  factors  for  converting  metric  (Sl)  units  of  measurement  to 
U.  S.  customary  units  and  U.  S.  customary  units  to  metric  (Sl)  units 
is  presented  on  page  3. 
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natural  coastal  configuration  where  water  is  free  to  move  artificially 
unobstructed  and  in  accord  with  the  forces. 

i* . The  use  of  any  of  the  previously  mentioned  open-coast  models 
in  simulating  circulation  and  water  level  conditions  is  made  difficult 
because  of  the  requirement  of  specifying  appropriate  boundary  conditions. 
This  is  especially  the  case  in  the  nearshore  region  where  these  models 
simulate  the  coast  boundary  as  an  infinitely  high,  continuous  wall. 

This  boundary  condition  neglects  overtopping  of  low-lying  land  and  bay 
communication  with  the  open  sea.  These  coastal  processes  have  a con- 
siderable effect  on  the  nearshore  water  levels  and  fluid  velocities. 
Furthermore,  it  is  precisely  this  zone  where  much  environmental  interest 
is  focused  and  the  use  of  numerical  models  should  prove  beneficial. 

5.  This  report  describes  a more  appropriate  open-coast  boundary 
condition  and  presents  results  from  11  verification  studies.  The 
boundary  condition  is  termed  the  finite  height  barrier  coast  (FHBC)  and 
is  incorporated  into  a two-dimensional  model  that  employs  an  orthogonal 
curvilinear  coordinate  system  with  telescoping  computing  cells.  This 
hydrodynamic  model  with  the  FHBC  is  named  SSURGE  III  and  uses  a slightly 

•5 

different  computing  scheme  than  that  previously  documented. 

6.  Similar  to  the  previous  model,  there  are  three  numerical 
programs  required  for  staging  the  production  of  surge  or  tide  computa- 
tions. The  first  numerical  program  determines  the  conformal  mapping 
coefficients  for  the  particular  study  region  of  the  continental  shelf. 

The  actual  coastline,  as  seen  on  National  Ocean  Survey  (NOS)  nautical 
charts,  is  smoothed  relative  to  portraying  small-scale  features.  This 
continuous  coastline  and  a deep  sea  boundary  curve  (perhaps,  following 
the  300-  or  600-ft  isobath)  are  the  input  to  the  conformal  mapping  pro- 
gram. The  coefficients  are  determined  and  inpu'  to  the  second  numerical 
program  which  determines  the  computing  grid  data  relative  to  a particular 
design  by  the  user.  The  grid  data  are,  in  part,  the  input  to  the 
SSURGE  III  program.  Other  input  includes  the  shelf  bathymetry,  FHBC 
data,  hydrograph  and  velocity  output  locations,  and  other  readily  de- 
terminable parameters.  In  actual  practice,  the  entire  process  is  not 
time-consuming  or  expensive. 


PART  II:  CONFORMAL  MAPPING 


7.  It  is  desired  to  conformally  map  a spatial  region  of  prototype 
space  of  the  continental  shelf  into  a rectangle  in  a mathematical  image 
plane,  in  which  the  coastline  and  deep  sea  boundary  curves  are  specifi- 
cally transformed  into  the  image  plane  as  constant  values  of  h Fig- 
ure 1 shows  the  details  of  the  transformation. 

3 

8.  It  can  be  shown  that  a conformal  transformation  satisfying 

the  above  constraints  is 

N 


x(C,ri)  = C 


(B  sinh  nkn  + C cosh  nkn)  sin  nkC 
n n 


y(C,n)  = B + 5 + y (B  cosh  nkn  + C sinh  nkn)  cos  nkC  (2) 

o n n 

n=l 

where 

k = u/X 

0£X;iX;0;iClX  (3) 

- 6 1 n e 

The  N values  of  B and  C and  B constitute  the  bi-curve  fitting 

n n o 

conformal  mapping  coefficients.  These  coefficients  are  determined  by 
matching  the  prototype  coastline  and  sea  boundary  curves  at  n = i 8 , 
respectively,  B also  being  a parameter  to  be  determined.  An  iterative 

3 

procedure  is  required  for  determining  6 and  the  coefficients. 
Essentially,  the  procedure  is  terminated  after  convergence  of  B and 
the  coefficients  or  after  the  transform-generated  coastline  and  sea 
boundary  curves  are  in  good  portrayal  of  those  specified  (personal  view 
of  the  user).  Typical  standard  deviation  between  transform-generated 
and  prototype  (slightly  smoothed  version)  over  the  length  of  the  curve 
is  of  the  order  of  2000  ft.  This  is  achieved  with  hO  iterations. 


For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  de 
fined  in  the  Notation  (Appendix  A). 


( 


Figure  1.  Conformal  transformation  mapping  prototype  space 

into  image  space 
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N = 100  and  at  a cost  of  9 min  of  CDC  6600  central  processing  time. 
Furthermore,  the  conformal  mapping  procedure  need  only  be  accomplished 
once  for  any  large  stretch  of  continental  shelf.  The  computing  grids 
appropriate  to  areas  within  that  region  are  generated  depending  on  the 
nature  of  the  particular  study,  using  the  mapping  coefficients  previously 
determined.  The  computer  cost  in  generating  the  grid  is  inconsequential 
to  the  personnel  time  required  in  the  grid  makeup  or  assigning  the  cell 
aver-aged  fluid  depth.  The  cost  in  generating  a grid  covering  a rela- 
tively small  extent  of  coast  must  be  weighed  against  a large  grid  of 
similar  high  resolution  over  a long  reach  of  coastline  which  costs  more 
per  prototype  hour  and  even  more  so  when  numerous  simulations  are 
required. 
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PART  III:  MATHEMATICAL  MODEL 


Stretched  Shelf  Coordinate  System 

9.  Consider  the  transform-generated  coastline  and  sea  boundary 
curves  shown  in  Figure  2.  The  orthogonal  curvilinear  mesh  associated 
with  the  (CjH)  coordinates  in  prototype  space  is  designated  the  shelf 
coordinate  system.  The  display  of  that  system  in  prototype  space  is 
shown  in  Figure  2a  and  its  image  in  Figure  2b.  Notice  that  only  a 
portion  of  the  entire  shelf  mapped  area  is  subsequently  employed  for 
long-wave  computations.  Furthermore,  the  calculations  cannot  readily  be 
performed  in  the  image  space  because  the  cells  are  unevenly  spaced. 
Centered  computations  are  a prerequisite  for  finite  differencing  of 
partial  differential  equations.  To  provide  an  evenly  spaced  computing 
grid  and,  at  the  same  time,  to  obtain  the  desired  spatial  resolution 
with  the  fewest  possible  computational  points  require  a second  trans- 
formation. This  transformation  preserves  the  orthogonal  property  and 
allows  for  the  independent  stretching  of  C and  of  n . The  grid  re- 
sulting from  the  second  transformation  is  the  evenly  spaced  computing 
grid.  Figure  2c.  The  coordinate  system  is  termed  the  stretched  shelf 
coordinate  system  (S,T). 

10.  The  stretched  shelf  coordinate  system  is  generated  by  in- 
dependently transforming  the  C and  n axes  in  the  following  manner: 

Given  the  nearshore  region  of  principal  interest,  the  values  of 
C along  the  coastline  are  determined  which  will  produce  a con- 
stant relatively  fine  increment  of  coastline  arc  length,  S 
In  this  area  of  prime  interest,  the  line  BC  in  Figure  2a,  ^ 

the  constant  increment  of  arc  length  is  equal  to  AS  . How- 
ever, regions  AB  and  CD  show  that  for  the  same  AS  as  above, 
there  is  a relative  expansion  of  the  increment  of  the  coastline 
arc  length.  The  functional  relation  between  5 and  S is 

S = S(S  (O)  ih) 

P 

where  the  expansion  of  Sp  with  respect  to  S is  specified 
by  an  (arbitrary)  expression  of  the  form 
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I 

i 

i 


s 

p 


'"a  + (S) 


(5) 


where  A , B,  and  C are  constants  for  each  expansion 
region  (k).  For  the  example  in  Figure  2,  there  are  three  ex- 
pansion regions.  The  ensemble  function,  k = 1,  2,  and 
3 , is  determined  such  that  at  interface  points  the  function 
is  continuous  and  has  continuous  first  derivatives. 

b.  Along  a particular  isoline  of  C » , the  values  of  n are 

determined  which  will  yield  a constant  change  in  the  time, 

AT  , required  for  a long  wave  traveling  at  the  local  free 
wave  celerity  to  proceed  from  the  sea  boundary  (-B)  to  the 

coast  (B).  The  long-wave  travel  time,  T , is  calculated 
■u  n ’ 

ty 


T 

n 


f ds 

S^(n) 


(6) 


where  is  the  distance  along  c'  , g is  the  acceleration 

due  to  gravity,  and  is  the  local  depth  relative  to  mean 

sea  level  for  a standard  basin.  This  procedure  provides  a 
fine  grid  spacing  in  the  nearshore  and  a coarse  grid  near 
the  deep  sea  boundary.  The  relation  between  n and  T 
is  given  by 


T = T(Tn(Sn(n)))  (?) 

The  incremental  values  of  T are  determined  from  Equation  6 
subject  to  the  (arbitrary)  expansion  relation  of  T(t  ) 
which  is  the  counterpart  of  Equation  5.  Actually,  the  rela- 
tion T(T  ) is  a convenience  (seldom  used  in  applica- 
tions) which  permits  an  additional  degree  of  freedom  in  ad- 
justing the  relative  spacing  between  isolines  of  n . In 
most  applications,  T = T and  the  value  of  AT  is  that 
which  divides  the  total  long-wave  travel  time  by  an  integer 
number  of  lines  of  n . The  selection  of  AT  relative  to 
the  coarse  deep  sea  spacing  is  based  upon  a compromise  for 
providing  adequate  resolution  of  the  hurricane  winds  and  the 
deepwater  surge  with  a minimimi  number  of  points.  At  first 
glance,  it  would  appear  that  the  T axis  for  hydrodynamic  long- 
wave calculations  is  time;  however,  this  is ‘not  the  case. 

Scale  factors  relating  the  transformation  of  n to  T are 
involved  in  the  equations  resulting  in  dimensions  of  length 
for  the  T independent  variab.'.e. 

11.  The  stretched  shelf  coordinate  system  provides  a grid  system 
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with  a finer  resolution  near  the  coast  than  at  the  deep  sea  boundary. 

The  expansion  curve,  3(3^),  stretches  the  alongshore  reach  of  the  grid 
while  maintaining  a finer  grid  in  the  area  of  principal  interest.  In 
this  manner,  an  economy  is  achieved  in  terms  of  the  number  of  grid  points 
required  for  long-wave  simulation.  However,  because  the  preferred  ex- 
pansion curves  dictate  the  locations  (in  prototype  space)  of  the  fluid 
depths  required  for  the  long-wave  calculations,  the  depth  field  must  be 
redefined  for  different  combinations  of  stretching  functions . 

Governing  Equations 

12.  The  vertically  integrated  form  of  the  quasi-linear  long-wave 
equations  in  a Cartesian  system  is  well  known. ^ The  appropriate  forms 
of  these  equations  in  the  stretched  shelf  coordinate  system  are 

Transport  (momentum)  equations 


^ - fO^  + It  (H  - H^)  = T,  - 


Fp  83 


8t  Fv  8T  ~ ’^T  ■ '’l 


Continuity  equation 


where  Q is  the  vol\ame  transport  per  unit  width  (units  of  length  /time), 
T is  the  wind  stress  divided  by  water  density,  a is  the  bottom 
frictional  resistance  stress  divided  by  water  density,  f is  the 
Coriolis  parameter,  g is  the  acceleration  due  to  gravity,  H is  the 
sea  surface  elevation  relative  to  mean,  sea  level,  D is  the  total 


12 


instantaneous  depth  of  water  (=H  - D^)» 


D the  water  depth  relative 


to  mean  sea  level,  and  is  the  hydrostatic  elevation  of  the  sea 

D 

surface  corresponding  to  the  atmospheric  pressure  anomaly.  The  in- 
dependent variables  are  time  (t)  and  the  coordinates  (S,T).  The  terms 
F , p and  v are  variable  scale  factors  associated  with  the  trans- 
formations. It  can  be  shown  that  F is  nondimens ional  and  given  by 

|l/2 

(11) 


F = 


The 


and  V scale  factors  are 


iL 

3S 


3S 

3S 


(12) 


and 


„ 3S  3T 
3n  n n 

^ “ 3S  3T  3T 
n n 


(13) 


It  can  be  shown  that  y(6S)  and  v(6T)  have  the  dimensions  of  length. 
2 

Furthermore,  F pv  corresponds  to  the  Jacobian  of  the  transformation 
in  the  sense  that 


R 


F pv  dS  dT 


(IM 


where  A is  the  area  of  a closed  region  in  prototype  space  whose 
corresponding  area  in  the  computing  space  is  R (see  Figure  2 for 
details).  This  relation  assures  that  all  area  enclosed  by  the  limits 
of  the  curvilinear  grid  in  prototype  space  is  accounted  for  in 
the  computing  space  by  use  of  the  scale  factors. 


13.  The  kinematic  wind-stress  components  Tg  and  are  re- 

to  their  x,y 
prototype  space  by 


lated  to  their  x,y  component  counterparts  (t  , t ) at  a given  point  in 

X y 
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T 


S 


T COS  6 + T sin  0 
X y 


(15) 


and 


T_  = - T sin  6 + T cos  0 
T X y 


(16) 


where 


(17) 


The  relation  between  the  wind  stress  and  wind  speed  at  a reference 
anemometer  level  (usually  taken  near  the  water  surface)  is  taken  as 


T 


(18) 


where  is  air  density  and  is  water  density,  is  a non- 

dimensioiial  drag  coefficient,  and  is  the  wind  speed  at  an  eleva- 

tion of  10  meters  above  the  water  surface.  The  value  of  K is  taken  as 


1.1  X 10 


-6 


, if  ^ 13.58  knots 


K = 


|i.l  + 2.5(1  - 13.58/W^Q^^lO  ^ , if  > 13.58  knots 


(19) 


* 3 5 ^ 

This  form  for  K was  used  in  previous  studies.  ’ Recent  compilation 

of  various  results  relating  to  the  wind  speed  at  10  meters  is  shown 

in  Figure  3.  As  reference,  the  dashed  line  in  Figure  3 shows  the 

relation  based  on  Equation  19  with  p /p  = 0.0012  . In  the  middle 

aw 

to  low  wind  speed  range  (say,  less  than  60  knots),  Equation  19  is 

representative  of  the  data.  It  is  this  range  of  wind  speeds  that  were 

used  by  Reid  and  Bodine^  in  studies  of  historical  storm  surges  in 

3 

Galveston  Bay.  Wanstrath,  on  the  other  hand,  in  studying  three  his- 
torical storm  surges,  each  in  different  coastal  regions  computed 
coastal  water  levels  in  good  agreement  with  observed  conditions  using 


SSURGEin 


WIND  SPEED.  KNOTS 


Figure  3.  Mean  values  of  the  neutral  drag  coefficient 
as  a function  of  low  and  high  wind  speeds  at  10  m height 
(from  Garratt,  1977)  and  that  used  in  SSURGE  III 


Equation  19  and  estimates  of  the  wind  furnished  hy  the  Hydrometeorologi- 
cal Section  (Hydromet)  of  the  National  Weather  Service  (NWS),  National 
Oceanic  and  Atmospheric  Administration  (NOAA).  Those  maximum  winds 
ranged  from  95  to  120  knots.  However,  the  current  understanding  of 
the  surface  winds  from  hurricanes  obtained  by  aerial  reconnaissance, 

anemometer  records  from  oil  platforms,  and  sophisticated  numerical 

T 8 

marine  boundary  layer  models  ’ is  that  the  surface  winds  at  10  meters 
are  not  as  severe  as  thought.  Consequently,  with  liberal  estimates  for 
the  maximum  winds,  a conservative  drag  coefficient  at  high  wind  speeds 
is  necessary.  The  procedure  does  provide  good  surge  results  as  evi- 
dencedby  these  and  previous  studies.  Actually,  all  previously  reported 
open-coast  surge  models  are  subject  to  recalibration  given  the  recent 
findings  of  hurricane  surface  winds. 

lU.  The  forms  of  the  seabed  frictional  resistance  terms  are 


and 


K Q 


'’t  2 S? 

D 


(20) 


(21) 


where 


(22) 


and  is  a variable  nondimensional  drag  coefficient  that  depends  on 
the  seabed  condition  and  water  depth.  For  typical  seabed  conditions, 

K ^ S(0..025,)  (23) 

° (1. 1*9)2 

where 

D*  = 0.25  , if  D < 0.5 

D*  = 1.0  , if  0.5  <.  D < 1.5  (2i*) 

D*  = 2.0  , if  1.5  1 D < 2.5 

and  all  depths  D and  are  in  feet. 
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Recursion  Equations 


15.  The  numerical  analogs  of  Equations  8-10  are  based  on  centered 

difference  approximations  of  all  terms.  The  algorithm  treats  the  time 
dependency  explicitly  and  employs  computing  lattices  as  shown  in  Fig- 
ure U.  The  recursion  equations  require  two  time  levels  to  complete  the 
cycle.  H values  as  well  as  and  x are  computed  on  the  lattice  at 

time  level  n . This  lattice  also  contains  the  permanent  storage  data 
of  D , 6 , X , and  Y where  X and  Y are  the  x and  y grid 
coordinates  of  H points.  For  a reference  index  (l,J),  the  transports 

Qg  and  are  computed  at  AS/2  and  AT/2  , respectively,  from 

the  H point  in  the  positive  axes  direction.  The  transports  are  com- 
puted on  the  lattice  at  time  level  n+1  . For  small  time  steps,  At  , 
it  is  inconsequential  that  x is  computed  at  n and  applied  at  n+1  . 
Furthermore,  x is  computed  at  H point  locations  and  averaged  with  its 
neighbor  in  the  S and  T directions,  respectively,  for  determination 
of  Xg  and  x^p  . These  approximations  for  x are  a result  of  efficient 
utilization  of  computer  time  and  memory  considerations  and  can  be  shown 
through  numerical  experiments  to  be  accurate. 

16.  Consider  that  the  transports  are  known  at  time  level  n-1 
and  the  H field  is  known  at  time  level  n-2  . The  recursion  formula 
for  interior  grid  H points  is 


H(l,J,n)  = H(l,J,n-2) 


2 At 

[Fg(I,J)Qg(l,J,n-l)  ■ 

,j)Qg(l-l,J,n-l)] 

Mp.ASpjp(l) 

m • 

Fj^(l,j) 

2 

2At 

[Fy(l,J)q.p(l,J,n-l)  ■ 

M 

1 

-1)Q^(I,J-I,n-1)] 

MpATVjj(j) 

5 — ' ■ a 

Fj^(l,j) 

2 

(25) 


where  1 , J , and  n Indices  express  the  S , T , and  time 
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1 


coordinates,  respectively,  At  is  the  numerical  time  step,  , 

and  F„  are  the  point  values  of  F from  Equation  11  determined  at 
n 

Qg  , and  H points,  respectively,  is  the  point  value  of  y from 

Equation  12  at  S^.  (typical  H grid  points  in  S direction),  v„  is 
i n 

the  point  value  of  v from  Equation  13  at  Tj  (typical  H grid  points 
in  T direction),  and  is  the  map  factor  relating  prototype  length 

(feet)  to  x,y  units.  The  simulation  is  initiated  at  n = 1 with  all 
transports  and  previous  H field  equal  to  zero. 

IT.  Given  the  H field  Just  computed  at  level  n and  the  previ- 
ous transports,  the  recursion  formulas  for  interior  grid  transports  are 


Qg(l,J,n+l) 


Qg(l,J,n-l)  2AtfQ.j, 


2AtgDg[H(l+l,J,n)  - H(l,J,n)  - Hg(l+l,J,n)  + Hj^(l,J,n)] 


M^ASGgyy(l)Fy(l,j)' 


and 


Q.p(l,J,n+l) 


+ Tg(l,J,n)  + Tg(l+1 

S L 


Q^(l,J,n-l)  2AtfQg 


.n)j 


(26) 


2AtgD^[H(l,J+l,n)  - H(l,J,n)  - Hg(l,J+l,n)  + Hg(l,J,n)] 


M^ATG^Vy(j)F^(I,J) 


At 


j^T5(l,J,n)  + Ti5(l,J+l,n)j 


(27) 


where 


Dg  = 0.5|H(l+l,J,n)  - D^(I+1,J)  + H(l,J,n)  - D^(l,J)J 


(28) 
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= 0.5[H(I,J+I,n)  - D^(I,J+1)  + H(l,J,n)  - (29) 

^ = 0.25|Qg(l,J+l,n-l)  + Qg(l,J,n-l) 

+ Qg(l-l,J-l,n-l)  + Qg(l-l,J,n-l)J  (30) 

^ = 0.25[Q^(l+l,J,n-l)  + Q^{l,J,n-l) 

+ Q^(l,J-l,n-l)  + Q^(l+l,J-l,n-l)J  (3l) 

“s  ■ 1 * (^)"|  ^^^1  (^sf  <32) 

“t  ' 3 + ■*  <33' 


, in  a similar  manner,  is  the  point  valUe  of  p from  Equation  12  at 
Sj  + AS/2  (typical  Qg  grid  points  in  the  S direction),  is  the 

point  value  of  v from  Equation  13  at  Tj  + AT/2  (typical  grid 

points  in  the  T direction),  tg  and  are  the  wind  stresses  in 

the  S , T directions,  respectively,  which  are  computed  at  H points, 

o rn 

and  are  the  seabed  drag  coefficients  computed  from  Equation  23 

with  D in  Equation  2h  given  by  Dg  and  , respectively. 

Boundary  Condition 

l8.  For  the  studies  reported  in  PART  IV,  two  general  types  of 
simulations  were  conducted.  Numerical  simulations  of  hydraulic 
physical  model  studies  require  that  the  lateral  boundaries,  1=1 
and  I = IM  , be  portrayed  as  walls.  The  sea  boundary  is  specified 
by  time-dependent  water  elevations,  H(l,l,n),  similar  to  that  observed 
in  the  physical  model.  In  some  of  these  studies,  reference  is  made 
to  hurricane  or  surge  simulations.  To  be  precise,  H(l,l,n)  is  a 
transient  response  function  but  no  wind  stress  is  applied  in  either  the 
physical  or  numerical  model.  The  coast  boundary  condition  is  the 
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Finite  Height  Barrier  Coast  (FHBC)  which  is  detailed  later. 

19.  Numerical  simulations  of  hurricane-induced  surges  require 
other  forms  for  the  lateral  boundary  conditions.  I^any  different  con- 
tinuum type  boundary  conditions  are  reported  in  the  literature.  In 
principle,  it  should  not  matter  which  of  these  is  used.  If  the  comput- 
ing grid  is  of  sufficiently  large  extent,  the  solution  in  the  central 
part  of  the  grid  near  the  coast  should  be  insensitive  to  the  particular 
one  used.  A popular  lateral  boundary  condition  which  is  also  used  in 
SSURGE  III  is 


(3M 


The  sea  boundary  condition  is  that  the  water-surface  elevation  is  placed 
in  equilibriiim  with  the  atmospheric  pressure  anomaly 

H(l,l,n)  = Hg(l,l,n)  (35) 

Along  the  lateral  and  sea  boundaries,  the  transports  are  computed  using 
modified  forms  of  Equations  30  and  31.  The  modification  is  to  compute 
an  average  value  using  only  the  two  adjacent  interior  points. 

20.  It  has  long  been  recognized  that  the  coast  boundary  condition 
plays  a major  role  in  the  nearshore  hydrodynamic  solution  of  the  govern- 
ing equations.  Most  open-coast  shelf  models  (whether  curvilinear, 
rectilinear,  finite  element,  implicit,  or  explicit)  treat  the  coast  as 
an  infinitely  high,  continuous  wall.  The  use  of  numerical  models  in 
simulating  circulation  and  water  level  conditions  is  made  difficult  be- 
cause of  application  of  such  a condition.  The  coast  boundary  condition 
developed  and  incorporated  into  SSURGE  III  is  the  FHBC . The  FHBC  pro- 
vides for  the  overtopping  of  low-lying  land  and  bay  communication  with 
the  open  sea. 

21.  The  FHBC  condition  is  to  allow  for  potential  ponding  areas 
landward  of  the  shoreline.  Figure  5 shows  the  conceptual  design  of  the 
barrier  coast  and  adjacent  bay.  The  bay  is  shown  to  be  of  several  incre- 
ments of  AS  in  length.  The  coast  flooding  routine  permits  a volme  of 
water  to  be  transported  across  the  nominal  shoreline  which  is  dependent 
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upon  barrier  heights,  predicted  water  level  at  the  coast,  channel  en- 
trance characteristics,  and  the  water  level  in  the  ponding  area, 

The  routine  is  applied  at  the  H time  levels  as  shown  in  Figure  h. 


22.  A prediction/correction  method  is  used  to  compute  H(l,JM,n)  . 
The  philosophy  is  taken  that  the  coast  barrier  as  viewed  by  the  long  wave 
is  essentially  a wall  over  much  of  the  time  of  surge  development,  re- 
flecting the  majority  of  the  long-wave  energy.  This  results  from  the 
consideration  that  long-wave  reflection  occurs  when  the  wave  experiences 

a significant  change  in  fluid  depth  over  distances  that  are  small  rela- 

9 10 

tive  to  the  wavelength.  ’ The  typical  nearshore  bathymetry  combined 
with  coastal  barrier  elevations  assures  that  the  coast  is  essentially 
a wall  which  leaks  water. 

23.  The  prediction/correction  method  is  to  first  assme  total 
reflection  at  the  boundary  and  then  to  correct  that  prediction  of  H 
based  on  a finite  height  coast.  In  this  manner,  the  method  provides  a 
smooth  transition  from  the  time-dependent  circumstances  where  the 
boundary  is  not  flooded  and  no  corrections  are  required  to  the  cata- 
strophic flooding  coast.  Other  relations  for  predicting  H at  the 
coast  than  that  of  total  reflection  were  considered.  In  particular, 

a water-surface  slope  projection  method  was  tested.  Results  of  these 
tests  which  are  presented  in  PART  IV  show  the  method  to  be  inap- 
propriate for  surge  simulations.  Other  tests  (results  not  presented) 
imply  that  the  slope  projection  method  is  applicable  to  tidal  (free 
wave)  simulations. 

2h.  At  time  level  n , the  water  level  at  the  coast,  h"*"  , is 
predicted.  If  the  water  levels,  h"*"  and  , exceed  the  minimum 

height  of  the  coast,  then  a volume  of  water  is  transported  across  the 
nominal  coast  for  the  time  interval,  n-1  to  n+1  , 

Submerged  Barrier 

'‘1  ■ 

where 

°h  = % 
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: 


(35) 

(36) 


(37) 


Cs  is  a nondimensional  submerged  barrier  coefficient  taken  as  O.U  , 

k k+1  k k+1 

and  are  lengths  of  the  coast  at  elevations  and  , 

respectively,  and  the  direction  of  the  flow  is  determined  by  the  sign 
of  (positive  means  water  is  removed  from  the  coastal  cell  and 

placed  in  the  ponding  area).  Barrier  heights  per  AS  section  centered 
on  H points  are  discretized  into  unit  elevations  (k).  The  length  of 
the  coast  at  each  elevation  is  obtained  from  topographic  maps,  beach 
surveys,  etc.  The  total  volume  of  water  crossing  the  submerged  coast 
is  the  summation  of  Equation  35  for  k = 0 , 1...K  where  K is  the 
integer  of  the  lesser  of  h"*"  or  • Considering  large  ponding  areas 
and  high  minimum  coastal  heights,  this  volume  is  generally  small 
compared  with  that  wnich  is  transported  through  channels  or  which  over- 
topped the  exposed  (on  one  side)  section  of  the  barrier.  These  submerged 
barrier  equations  reflect  the  assumption  that  frictional  effects  are 
dominant  in  the  bottom  layer  which  is  continuous  across  the  boundary. 

25.  For  the  overtopped  barrier,  the  volume  of  water  crossing 
the  coast  is 

Exposed  (on  one  side)  Barrier 


where 


- , If 

- 2^  . if  Ht  < i 


(38) 


(39) 


C„  is  a nondimensional  exposed  barrier  coefficient  taken  as  0.2  ai.d  the 

rj 

direction  of  flow  is  determined  by  the  sign  of  . The  total  volume 

of  water  overtopping  the  coast  is  the  summation  of  Equation  38  for 

+ 

k = 0,1, ...M  where  M is  the  integer  of  the  greater  of  H or  . 


2h 


26.  The  volume  of  water  associated  with  channel  communication  is 


= + 2At 


- "kl  * ''c”oV^c) 


>S,I  • if 

, if 

, if 
, if 


h"^  > 0 
> 0 

h"^  < 0 
< 0 


and 

and 

and 

and 


\ > 0 
\ < 0 

«b  " ° 


(1*0) 


(hi) 


is  the  channel  width,  is  a nondimensional  channel  discharge 

coefficient  determined  for  each  entrance,  A is  the  channel  cross- 
sectional  area  at  mean  sea  level,  and  the  sense  of  the  flow  is  taken 
toward  the  low  head  side.  An  estimate  of  CI^A^  is  possible  if  tidal 
observations  are  available  inside  and  outside  a bay  with  a constricted 
opening  to  the  sea 


A /32ttH 

C A = 2. 

DC  3gBp 


(U2) 


where  A^  is  the  surface  area  of  the  bay  at  mean  sea  level,  T is 
the  tidal  period,  is  the  tidal  range  outside  the  bay,  and  is 

a nondimensional  bay  parameter  that  is  dependent  on  the  response 

or  phase  lag.  H.  is  the  tidal  range  inside  the  bay.  The  bay  parameter 

^ 11 
as  a function  of  response  or  phase  lag  is  shown  in  Figure  6. 

27,  Consider  that  h"*"  values  along  the  coast  (uncorrected  for 

flooding)  are  known  at  time  level  n and  and  A^  are  known  at 

level  n-2  where  A^  is  the  bay  storage  area.  The  volinne  of  water 

entering  (or  leaving)  the  ponding  area  is  determined  for  each  segment 

k k 

by  the  appropriate  sum  of  , and  . The  ratio  of  that 

volume  to  the  surface  area  of  the  grid  block  representative  of  H , 

Aj  , provides  the  incremental  correction  to  the  predicted  coastal  water 
level  to  conserve  mass.  This  value  is  stored  in  H(l,JM,n)  for  use  in 
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RESPONSE  FACTOR,  H /H 


BAY  parameter, 


* 

NOTE:  H i = TIDAL  RANGE 

INSIDE  THE  BAY 

H = TIDAL  RANGE 

OUTSIDE  THE  BAY 

Figure  6.  Bay  parameter  Bp  as  a function  of  response  factor 

or  phase  lag 
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PHASE  LAG,  DEG 


the  transport  computations  at  level  n+I  . The  new  value  of  at 

time  level  n is  at  n-2  plus  the  incremental,  water  level  change 

from  the  instantaneous  and  even  distribution  of  the  entering  volume 
(the  s\im  for  all  appropriate  AS  segments)  on  at  level  n-2  . 

Prom  at  n , a new  bay  storage  area  is  obtained.  The  above  is 

repeated  for  each  ponding  area  and  its  designated  coastal  segments. 

Hurricane  Wind  and  Atmospheric  Pressure  Models 

28.  The  hurricane  wind  and  atmospheric  pressure  models  as  given 
12 

by  Jelesnianski  are  employed  in  the  surge  studies  with  some  minor  ad- 
ditions. The  x,y-wind  components  for  a stationary  storm  at  the  H grid 
points  are 

\ = [-(Xg  ••  Xg)  sin  (t>  - (Yg  - Y^)  cos  d]F(rjj)  (1*3) 

and 

w = — ["(X  - X ) cos  (|)  - (Y  - Y ) sin  (^1p(r„) 

Se  JH 


where 

r„-  ttM 


(1*5) 


W is  the  stationary  storm  maximum  wind,  <(i  is  the  wind  ingress  angle 
n 

reflecting  the  inward  flow  relative  to  that  wind  vector  tangent  to  the 

isovel,  R„  is  the  distance  from  the  storm  center  (X  ,Y  ) to  region  of 
H 6 6 

maximum  winds,  and  X and  Y are  the  (x,y)  coordinates  of  H grid 

g g 
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points.  The  translation  of  the  storm  provides  an  alteration  in  the  wind 
field.  The  x,y  supplemental  wind  components  due  to  storm  translation 
are 

^ (h6) 

XX  y n 

and 

t = (-T  sin  + T cos  i())G{r  ) (U7) 

y X y H 

where 

i|;  = a - <t>  - 90  (I48) 


if 


if 


r > R 
H - H 


(U9) 


T and  T are  x,y  components  of  the  forward  speed  of  the  storm  and 
X y 

a is  a rotation  angle  used,  primarily,  to  fine  tune  the  alignment  of 
maximum  winds.  For  standard  operation  of  the  wind  model,  a is  set 
equal  to  90°  + (p  . The  x,y-wind  stress  components  for  a moving  storm 
are 


and 


where 
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and  K is  given  by  Equation  19-  The  stress  components  in  the  stretched 
shelf  coordinate  system  are  determined  by  applying  Equations  50  and  51 
in  Equations  15  and  l6.  In  an  alternate  manner. 


W„  = W cos  0 + W sin  0 
S X y 


(5M 


and 


Wn,  = -W  sin  0 + W cos  0 
T X y 


(55) 


where  0 is  given  by  Equation  IT  and  Wg  and  are  the  wind  compo- 

nents in  the  (S,T)  system.  The  stress  components  are 


and 


^ 


Tt  = k 


/ 2 2 \ 
(Wg  . 


(56) 


1/2 


(5T) 


29.  The  surface  atmospheric  pressure  field  associated  (but  not 
dynamically  coupled)  with  the  hurricane  wind  is 


p = p + (p  - p )e 
o ® o 


-R„/r 


(58) 


where  is  the  central  pressure  and  is  the  far-field  pressure. 

It  is  often  observed  that  relatively  high  winds  remain  along  the  coast 

after  the  storm  has  proceeded  inland.  Through  numerical  experiments, 

it  is  found  that  reasonably  good  comparison  between  observed  (Hydromet) 

and  computed  winds  is  obtained  by  setting  to  be  the  distance  the 

n 

storm  center  is  from  the  coast,  specifying  Wp  to  be  desired  winds  at 
the  coast  and  setting  a (measured  clockwise  from  the  storm  movement) 
such  that  the  maximum  wind  region  is  along  the  coast.  This  procedure 
allows  slightly  longer  simulations  to  be  performed  after  landfall  than 
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would  be  permitted  otherwise  due  to  poor  wind  field  portrayal.  However, 

R„  in  the  pressure  expression  (Equation  58)  is  not  allowed  to  increase, 
n 

If  this  procedure  is  followed  in  the  numerical  program  (it  is  optional), 
Rjj  is  held  constant  at  a value  RHIT  in  nautical  miles  (n.mi.)  after  the 
time  of  storm  landfall  (THIT)  which  is  input  in  hours  after  the  start  of 
the  simulation. 

30.  The  various  program  parameters  that  are  used  to  generate  the 

13 

wind  and  atmospheric  pressure  fields  for  Hurricanes  Donna,  Flossy, 

lU  15 

Hilda,  Betsy,  and  Carmen  are  given  in  Tables  1-5.  Only  for  the 
postlandfall  winds  from  Hurricanes  Hilda,  Betsy,  and  Carmen  were  the 
normal  routines  modified  in  the  above  manner.  This  is  seen  by  noting 
the  number  for  THIT. 

31.  The  term  H^  in  Equations  8,  9,  26,  27,  and  35  is  computed 

D 

in  feet  of  water  by 

Hg  = 0.0328(P^  - P)  (59) 

where  P , named  PINF  in  Tables  1-5,  and  P are  in  pressures  of 

00  ' O 

millibars. 

32.  For  historical  storms,  all  necessary  input  parameters  can  be 

determined  if  Hydromet  has  assessed  the  maximum  winds  following 

the  same  procedures.  This  procedure  probably  results  in  a liberal  esti- 
mate for  . Storm  surge  simulations  for  forecasting  purposes  typi- 

K 

cally  involve  storm  parameters  of  track,  P„  , , and  R^  . The 

only  other  necessary  input  for  this  wind  model  is  Wg  . This  can  be 
determined  from  the  correlation  of  a few  severe  historical  storms  where 
Hydromet  has  estimated  W and  that  value  which  is  predicted  by  the 
Standard  Project  Hurricane  (SPH)  method  using  only  the  required  fore- 
cast storm  parameters.  A cursory  examination  shows  that  the  ratio  Wg 
to  probably  in  the  range  1.2  to  1.3.  This  approach  is 

suggested  only  as  a temporary  solution.  Actually,  a dynamic  marine 
boundai'y  layer  hurricane  model  is  needed  that  is  independently  cali- 
brated to  several  storms  of  record  and  verified  against  even  more 
historical  storms.  This  wind  model  is  comparable  to  the  treatment  of 
the  physics  of  the  atmosphere  as  SSURGE  III  is  to  the  hydrodynamics. 
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PART  IV:  MATHEMATICAL  MODEL  APPLICATIONS 


1 

I 

i 

I 

33.  The  results  of  five  studies  are  reported.  In  all  but  one 
study,  simulations  of  different  events  were  performed  to  assess  the  be- 
havior of  the  FHBC  under  varying  conditions.  Two  studies,  the  i 

Galveston  Bay  and  Murrells  Inlet  studies,  compare  the  free  wave  solution  j 

as  observed  in  the  physical  model  and  computed  in  the  mathematical  model.  j 

Furthermore,  the  solutions  are  for  the  nearshore,  open  sea  waters  and  i 

not,  per  se,  for  the  interior  of  the  bay  or  inlet.  In  some  cases,  the 
prototype  data  used  to  calibrate  the  physical  model  are  used.  The 
plates,  reflecting  the  results  of  these  studies,  refer  to  "surge"  or 
"hurricane"  simulations.  As  noted  previously,  the  sea  boundary  water 
level  is  specified  as  a transient  response  function  and  no  wind  stress 

is  applied  in  either  the  physical  or  numerical  model.  ■ 

3^.  The  remaining  three  studies  deal  with  a computing  grid  sensi- 
tivity study  and  actual  hurricane  surge  simulations  from  historical 
storms.  The  grid  sensitivity  study  compares  the  surge  results  obtained 
from  variable  spaced  and  evenly  spaced  rectilinear  computing  grids . The 
storm  surge  from  Hurricane  Donna  is  simulated  along  the  west  coast  of 
Florida  for  a wall  and  flooding  coast  conditions.  The  final  study  is  of 
the  storm  surges  from  Hurricanes  Flossy,  Hilda,  Betsy,  and  Carmen  which 
affected  the  Louisiana  coast  from  Atchafalaya  Bay  to  the  Mississippi 
River. 

Galveston  Bay  Physical  Model  Study 

IT 

35.  The  Galveston  Bay  physical  model  was  constructed  for  study 
of  hurricane  surge  routing.  The  physical  model  results  or  observed 
(prototype)  data  are  from  that  report.  The  general  computing  grid, 
location  map  of  tide  or  velocity  gages,  particular  computing  grid  points, 
and  channel  entrance  details  are  presented  in  Plate  1.  The  FHBC  data 
are  given  in  Table  6. 

36.  The  procedures  used  in  both  the  numerical  and  physical  models 
were  to  input  a water  level  at  the  respective  sea  boundaries  such  that 
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the  hydrograph  at  Pleasure  Pier,  sta  1,  is  reproduced.  This  pro- 
cedure is  followed  for  the  astronomical  tide.  Hurricane  Carla  and  large 
radius,  high  translation  (LRHT)  hurricane  simulations.  The  hydrograph 
at  Pleasure  Pier  from  the  numerical  model  and  the  prototype  (essen- 
tially, that  recorded  in  the  physical  model)  is  presented  in  Plate  2.  i 

Having  achieved  the  proper  input,  the  numerical  portrayals  for  the  flow  j 

in  the  nearshore,  channel  entrance  velocity,  and  hydrograph  are  \ 

obtained.  Plates  3 and  U show  the  depth-averaged  velocity  fields,  as 
determined  by  SSURGE  III,  in  the  Galveston  Channel  entrance  region  for 
the  flood  and  ebb  stages,  respectively.  For  the  tide  simulation,  the 
numerical  program  was  specially  altered  to  crudely  portray  the  north  and 

south  entrance  Jetties.  By  doing  this,  better  agreement  in  the  compari-  i 

son  between  the  prototype  and  numerical  model  entrance  channel  velocity 
is  accomplished.  This  velocity  comparison  is  shown  in  Plate  5 where,  in 
addition,  the  velocity  at  Rollover  Pass  is  present.  For  the  tide  range. 

Rollover  Pass  is  essentially  a wall.  Consequently,  one  can  conclude 
that  if  a wall  were  placed  across  the  main  entrance — not  following  the 
FHBC — an  inappropriate  simulation  would  result.  The  response  of 
Galveston  Bay  at  selected  "mean"  stations  and  the  computed  H^  are 
shown  in  Plate  6.  Although  the  velocity  fields  shown  in  Plates  3 and  U 
are  seemingly  very  representative  of  what  one  should  expect  to  occur, 
the  detailed  entrance  velocity  comparison  shows  that  the  FHBC  overreacts 
to  the  small  semidiurnal  component.  This  should  be  expected  in  such 
cases  because  the  numerical  solution  does  not  contain  the  inertia  of 
Galveston  Bay.  Considering  that  the  computer  cost  for  a simulation  is 
about  a couple  of  dollars,  reasonable  results  are  obtained. 

3T.  The  hydrographs  at  Pleasure  Pier  for  the  Hurricane  Carla 
simulation  are  shown  in  Plate  T.  The  FHBC  solution  at  the  grid  boundary 
for  the  channel  (Plate  8)  shows  a small  O.U-ft  drawdown  relative  to  the 
peak  of  sta  1.  The  computed  peak  elevation  occurs  nearly  in  phase  with 
sta  1.  For  comparison,  the  hydrograph  at  sta  h is  presented  and  shows 
it  to  lag  that  of  the  open-coast  records.  Consequently,  it  is  thought 
that  the  numerical  solution  at  the  boundary  does  follow  that  expected 
for  the  nearshore  open-coast  waters.  No  physical  model  entrance 
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channel  velocities  were  recorded.  The  numerical  model  velocities  of  the 
Galveston  Entrance  Channel  and  Rollover  Pass  are  presented  in  Plate  9 
for  comparative  purposes  with  the  tidal  results  and  the  LRHT  hurricane 
results.  The  computed  is  in  reasonable  agreement  with  the  bay 

hydrograph  at  sta  T shown  in  Plate  10. 

38 . The  special  procedure  for  numerically  simulating  the  jetties 
was  not  included  in  the  program  for  this  or  the  LRHT  case  because  the 
jetties  were  readily  overtopped  at  these  water  levels.  It  is  not  known 
if  the  channel  velocity  should  be  amplified  in  accord  with  that  found 
for  the  tidal  simulation.  It  is  recommended  that  the  numerical  program 
be  generalized  such  that  it  can  treat  open-coast  structures  of  finite 
height.  In  the  study  of  Murrells  Inlet,  it  is  difficult  to  assess  the 
model's  behavior  because  of  the  lack  of  such  program  features. 

39.  Plates  11-lL  are  for  the  LRHT  simulation  and  follow  in  similar 

designation  those  for  the  Hurricane  Carla  simulation.  Again,  it  is  noted 
that  the  FHBC  solution  at  the  grid  boundary  shows  a drawdown  of  O.T  ft 
relative  to  the  peak  at  sta  1.  The  phase  of  the  peak  occurs  nearly 
midway  between  sta  1 and  U (Plate  12).  However,  relative  to  the 
amplitudes,  the  FHBC  boundary  solution  would  appear  to  follow  that  of 
the  expected  entrance  condition.  The  comparison  of  Plates  9 and  13  pro- 
vides qualitative  assurance  that  the  model  is  responding  in  accord  with 
expected  conditions.  Although  the  agreement  between  and  sta  7 in 

phasing  and  amplitude  is  not  as  good  as  in  the  other  simulations,  H^ 

is  representative  of  the  bay  water  elevations  (Plate  lU). 

Murrells  Inlet  Physical  Model  Study 

1 8 

Uo.  The  physical  model  study  of  Murrells  Inlet  was  recently 
completed.  The  physical  model  results  for  the  astronomical  tide  simula- 
tion are  from  that  report.  The  same  entrance  details  (Plan  IH)  were 
used  for  the  tide  and  physical  model  surge  simulations.  The  two  surge 
experiments  were  conducted  by  this  author  with  assistance  from  others 
in  the  Hydraulics  Laboratory. 

1*1 . The  computing  grid  and  location  map  of  tide  gages  are 
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presented  in  Plate  15.  The  FHBC  data  are  shown  in  Table  7.  The  physical 
model  entrance  details  and  the  computing  grid  for  that  area  are  shown  in 
Plate  l6.  The  numerical  program  was  not  specially  altered  to  reflect 
these  finite  height  open  sea  structures  (jetties,  weirs,  or  training 
dikes).  As  noted  previously,  the  program  can  readily  be  generalized  for 
proper  treatment  of  such  situations.  The  results  comparing  the  physical 
and  numerical  models  are  presented  here  for  two  reasons.  First,  Murrells 
Inlet  offers  a contrasting  study  to  the  much  larger  Galveston  Bay; 
second,  data  will  exist  for  a subsequent  study  with  an  improved  model. 

h2.  Tide  sta  10  provides  the  sea  boundary  forcing  condition  to 
the  mathematical  model  (Plate  IT).  Essentially,  it  is  also  the  control 
hydrograph  for  the  physical  model  where  water  is  added  or  removed  at  the 
head  bay  (approximately,  Y = -10  , Plate  15)  such  that  the  control 
hydrograph  is  followed.  The  flood  and  ebb  stages  of  the  depth-averaged 
velocity  fields  for  the  numerical  tide  simulation  are  shown  in  Plates  l8 
and  19,  respectively.  The  velocity  fields  are  illustrative  and  repre- 
sentative for  the  portrayed  situation.  The  time  history  of  the  numerical 
entrance  velocity  located  in  close  proximity  to  velocity  sta  1 
(Plate  l6)  is  compared  in  Plate  20  with  that  recorded  from  the  physical 
model.  The  flood  velocity  from  the  nimierical  simulation  agrees,  probably 
fortuitously,  with  the  physical  model  results;  however,  the  ebb  veloc- 
ity comparison  does  not.  It  is  noted  that  these  results  for  the  flood 
and  ebb  velocity  comparison  also  follow  for  the  two  surge  simulations. 
Plate  21  compares  the  observed  interior  hydrographs  and  computed  . 

Clearly,  a representative  time-dependent  transport  of  water  into  the 
bay  is  achieved  through  the  FHBC  procedure. 

^3.  Plates  22-25  compare  the  results  for  a 9- ft  surge  simulation. 
The  sea  boundary  hydrograph  is  presented  in  Plate  22  where  the  surge 
height  is  actually  9.8  ft  relative  to  mean  low  water  (mlw).  The  name, 
9-ft  surge  simulation,  is  carried  for  convenience.  The  comparison  of 
the  momerical  and  physical  models'  water  levels  at  the  entrance  is 
shown  in  Plate  23.  The  drawdown  through  the  entrance  channel  is 
represented  in  the  results  of  the  numerical  model.  The  velocity  com- 
parison (Plate  2h)  shows  fair  agreement  in  the  flood  and  poor  agreement 
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in  the  ebh.  Plate  25  shows  that  the  computed  peak  bay  water  level  is  in 
agreement  with  several  peak  water  levels  at  selected  stations;  but 
decreases  too  rapidly  relative  to  observed  conditions. 

I+I4.  Plates  26-29  show  the  results  for  12- ft  surge  experiments. 

The  coastal  and  inland  topographies  of  the  physical  model  were  accu- 
rately portrayed  to  an  elevation  of  +10  ft  (relative  to  mlw).  Flooding 
in  excess  of  10  ft  is  constrained  by  the  walls  of  the  physical  model. 

U5.  The  sea  boundary  hydrograph  (Plate  26)  shows  the  results  of 
a n\americal  experiment  which  in  addition  to  specifying  the  (known) 
forcing  condition  (H^q)  , allows  that  input  to  be  instantaneously  ad- 
justed for  the  outward  transport  of  water  at  the  boundary.  This  numeri- 
cal boundary  condition  is  in  close  analog  with  that  of  the  physical 
model  controller  which  adds  or  removes  water  given  the  control  hydro- 
graph and  that  continuously  sensed  at  the  gage.  The  sea  boundary 
condition  is 


Hio(n) 


Hjo(n)  » 


‘'gcT 


if  Qrp  ^ 0 (inward  flow) 


if  Qrp  < 0 (outward  flow) 


(60) 


where  is  the  observed  hydrograph  at  sta  10  and 


,1/2 


(^gDg  j is 

a constant  wave  speed  appropriate  for  deep  water.  For  this  application, 
D is  taken  as  500  ft.  The  second  term  in  Equation  60  represents 

D 

an  additional  rise  in  the  water  level  from  the  outward  propagating  wave 
energy.  The  homogeneous  form  of  Equation  60  is  most  pleasing  since 
it  involves  both  the  water  level  and  the  local  computed  transport.  In 
fact,  without  applying  this  condition  the  computations  were  unstable 
in  the  decreasing  stage  of  the  surge  where  the  time  step  was  taken  at 
10  sec.  Pseudostable  computations  resulted  when  the  time  step  was 
reduced  to  6 sec  as  shown  by  the  dashed  line  in  Plate  27.  The  solid 
line  in  Plate  27  shows  the  results  with  Equation  60  and  the  larger  time 
step  of  10  sec.  It  is  seen  that  without  the  radiation- type  sea 
boundary  condition,  significant  trapping  of  wave  energy  occurs  between 


1 
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the  coast  and  sea  boundary.  The  energy  directed  toward  the  sea  boundary 
in  the  decreasing  stage  is  much  greater  in  this  experiment  than  in  the 
9- ft  surge  experiment  because  of  the  overtopping  of  the  entire  coastline. 
For  this  reason,  there  was  found  to  be  no  practical  difference  in  the 
9- ft  surge  results  using  Equation  60  and  those  previously  presented. 

More  testing  of  the  radiation  condition  is  required  before  it  can  be 
accepted  for  general  use. 

U6.  Similar  results  concerning  the  drawdown  through  the  entrance 
channel  as  shown  in  Plate  23  are  obtained  in  the  12-ft  surge  simulations. 
The  previous  discussions  comparing  the  entrance  channel  velocities  and 
inland  hydrographs  for  the  9-ft  surge  experiments  are  applicable  also  in 
this  case.  Plate  28,  which  compares  the  entrance  channel  velocities  for 
the  physical  and  numerical  models,  shows  fair  agreement  in  the  flood 
stage  and  poor  agreement  in  the  ebb  stage.  Plate  29  shows  the  computed 
and  selected  inland  hydrograph  stations.  As  before,  decreases 

too  rapidly  relative  to  observed  conditions. 


i+T.  Two  hypothetical  storms.  Hurricanes  A and  B,  are  employed  to 
induce  surges  in  a standard  (or  fixed)  basin  of  specified  length,  width, 
and  bathymetry.  Hurrican  A is  an  onshore  moving  storm  while  the  track 
of  Hurricane  B is  parallel  to  the  coast.  A varic-ty  of  rectilinear  com- 
puting grids  are  used  to  portray  the  standard  basin.  Grids  A,  B,  and  C 
are  unevenly  spaced  in  x and  y,  but  they  all  have  the  same  nimiber  and 
alignment  of  points  in  the  x-direction.  Each  of  these  grids  have 
different  numbers  of  points  in  the  y-direction  since  the  same  basin 
width  is  used.  Grid  D is  an  evenly  spaced  grid  in  both  the  x-  and 
y-directions . Plate  30  shows  the  standard  basin,  the  various  grids, 
the  storm  parameters,  and  hurricane  tracks.  In  particular,  grids  A,  B, 
and  C have  identical  computing  points  at  y = 82  n.mi  (or,  8 n.mi.  from 
the  coast). 

I48.  There  are  no  analytical  or  known  solutions  such  that  an  as- 
sessment can  be  performed  of  the  accuracy  of  the  results  from  one  grid 
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relative  to  the  results  from  another  grid.  Rather,  the  approach  taken 
here  is  to  show  that  the  same  results  are  obtained  from  the  variable 
spaced  grids  which  involve  the  stretching  functions  (ij,v)  and  the  evenly 
spaced  grid.  Results  from  grids  B,  C,  and  D are  presented  in  Plate  31 
for  an  onshore  storm  and  in  Plate  32  for  an  alongshore  storm.  The 
results  from  the  unevenly  spaced  grids  for  each  type  track  are  equiva- 
lent to  that  from  the  evenly  spaced  grid.  The  equivalent  computations 
from  grids  B and  C should  not  be  generalized  to  typical  open-coast 
surge  conditions  where  varying  bathymetry  and  coastal  configurations 
will  play  a role  in  the  calculations.  The  term  "coastal  surge  en- 
velope" refers  to  the  peak  surge  water  elevations  along  the  coast  with- 
out regard  to  the  time  of  its  occurrence.  All  simulations  are  for  a 
22-hr  period  of  time  with  landfall  for  Hurricane  A at  l8  hr.  Hurricane 
B,  on  the  other  hand,  is  locatea  only  30  n.mi.  from  the  right  lateral 
boundary  where,  at  the  start  of  computations,  hurricane  force  winds  are 
imparted  to  the  basin  (initial  basin  water  conditions  are  H , Qg  , 
and  Qrp  = O).  This  procedure  is  not  recommended  for  general  surge 
modeling  since  the  basin  water  is  not  permitted  an  initial  transition 
period  for  spin-up.  Consequently,  transient  conditions  cause  the  surge 
envelope  bubble  exhibited  in  Plate  32  for  all  grids  at  x = 3^0  n.mi. 

The  envelope  is  representative  of  that  expected  (thau  is,  a constant) 
in  the  range  x = 120  n.mi.  to  x = 270  n.mi.  This  stretch  of  coastal 
surge  development  corresponds  to  the  last  10  hr  of  the  surge  simulation. 
It  may  be  possible  to  minimize  the  initial  model  spin-up  period  and,  at 
the  same  time,  reduce  the  ti-ansient  bubble.  Development  of  such  pro- 
cedures would  result  in  significant  computer  cost  savings. 

kg.  For  the  linear,  free  wave  simulation  in  an  enclosed,  constant 
depth  body  of  water,  it  can  be  shown  that  the  niunerical  analogs  of  the 
governing  equations  of  motion  conserve  mass  and  energy.  As  noted  in 
PART  III,  the  coastal  water  level  is  first  predicted  (and  subsequently 
corrected  for  flooding)  by  assuming  total  energy  reflection  at  the 
boundary.  This  condition  is  also  used  in  the  above-mentioned  free  wave 
seiche  simulations.  An  alternate  coastal  water  level  prediction  method 
is  studied  which  linearly  pro.lects  the  water  level  to  the  coast  boundary 
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from  the  slope  of  the  water  surface.  This  boundary  condition  does  not 
guarantee  conservation  of  mass.  Comparison  of  surge  results  from 
grids  A,  B,  and  C using  the  slope  projection  boundary  condition  and 
grid  B using  the  mass  continuity  (total  energy  reflection)  method  are 
presented  in  Plates  33  and  3^  for  Hurricanes  A and  B,  respectively. 

Even  for  grid  A,  which  has  a very  small  spacing  in  the  y-direction  at 
the  coast,  the  slope  method  proves  inappropriate  for  surge  computations . 
Not  only  are  poor  results  computed  at  the  coast  in  the  principal  surge 
development  area  but,  since  the  equations  of  motion  are  integrated  in 
time,  poor  agreement  in  the  hydrographs  is  observed  at  8 n.mi.  from  the 
coast  as  shown  in  Plate  35.  However,  in  other  applications  the  slope 
projection  boundary  condition  did  yield  results  (not  presented)  similar 
to  those  previously  shown  for  the  Galveston  Bay  studies.  Qualitatively, 
reasonable  results  are  seen  to  occur  when  the  flooding  coast  is  used 
with  the  slope  projection  method  in  grid  B for  the  onshore  storm 
(Plate  36)  and  the  alongshore  hurricane  (Plate  37).  Further  testing  of 
projection  method  is  warranted  centering  on  nonlinear,  quadratic,  or 
cubic  projections.  The  method  may  prove  useful  in  other  applications 
where  projection  to  a boundary  is  required. 


West  Coast  of  Florida,  Hurricane  Donna  Surge  Study 

50.  Hurricane  Donna,  which  occurred  in  September  I96O,  paral- 
leled the  southwest  coast  of  Florida,  crossing  the  Florida  Keys  and 
landfalling  near  Cape  Sable.  Table  1 presents  the  storm  parameters. 
Plate  38  shows  the  location  map  of  the  storm  track  for  the  last  I8  hr 
of  the  30-hr  simulation  and  the  computing  grid.  The  Florida  Keys  are 
portrayed  as  a wall  (Plate  38)  because  of  the  land  and  the  highway  ele- 
vations of  U.S.  1.  From  the  surge  results  and  observed  conditions, 
portions  of  the  highway  in  the  central  keys  were  washed  out  or  over- 
topped; but  no  attempt  is  made  to  account  for  this. 

51.  The  FHBC  condition  is  applied  to  the  coast  from  10  n.mi. 
north  of  Cape  Romano  to  the  end  of  the  grid  in  i’lorida  Bay.  Table  8 
presents  the  FHBC  data.  Particular  to  this  coastal  region  is  the 
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Everglades  National  Park  and  the  dense  population  of  both  red  and  black 

mangroves.  These  trees,  essentially,  are  the  coastal  barrier  rather 
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than  sand  dunes,  sea  walls  etc.  They  are  described  in  many  studies 
and  were  described  through  personal  communications  with  Mr.  Ben 
McPherson,  Aquatic  Biologist,  U.  S.  Department  of  Interior,  Geologic 
Survey,  Water  Resources  Division,  Miami,  Fla.  The  typical  widths 
of  the  mangrove  stand  are  from  1 to  10  miles  and  the  minimum  heights 
range  from  10  to  15  ft,  increasing  to  25  ft  in  more  inland  areas.  The 
red  mangrove  has  a thick  tangle  of  prop  roots  which  act  to  catch  and 
hold  sediment.  An  agile,  6- ft  tall  person  can  supposedly  move  about  on 
these  roots,  traversing  the  mangrove  stand  from  the  coast  inland  to  the 
large  relatively  flat  swamp  and  prairie.  Above  this  is  a thick,  inter- 
weaving canopy  of  tree  limbs,  branches,  and  leaves.  Since  the  mangrove 
characteristics  relative  to  long-wave  attenuation  are  not  known,  though 

it  would  appear  that  little  or  no  attentuation  should  exist  for  the 
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trunk  portion  of  the  tree  stand,  the  approach  taken  in  this  study  is 
to  assume  a standard  tree  barrier  coast.  This  tree  coast  is  a barrier 
from  msl  to  3 ft  (a  nominal  beach  berm  plus  the  prop  root  system),  from 
3 to  9 ft  is  the  trunk  portion  where  no  attenuation  is  assumed,  and 
from  9 to  12  ft  is  another  barrier  portraying  the  canopy.  Only  a small 
length  of  coastline  on  Cape  Sable  is  free  of  the  mangrove  as  seen  in 
Table  8.  The  normal  FHBC  routine,  including  barrier  coefficients,  was 
followed  except  in  the  situation  where  if  the  predicted  coastal  water 
level  opposite  a tree  barrier  exceeds  9 ft,  the  water  level  was  taken 
as  9 ft  for  flooding  computations  appropriate  to  the  trunk  portion. 

For  example,  if  the  predicted  coastal  water  level  is  13.5  Tt,  only  9 ft 
is  used  for  flooding  the  trunk  portion  of  the  tree  barrier  and  1.5  ft 
is  used  for  overtopping  the  canopy. 

52.  There  are  four  discrete  ponding  areas  whose  storage  area  in- 
formation is  presented  in  Table  8.  Flood  Region  I is  from  Naples, 

Fla.,  to  Marco  Island  and  bounded  inland  by  highways  U.S.  *4l  and  Fla.  92. 
Flood  Region  II  is  from  Marco  Island  to  Everglades  City,  Fla.,  and 
bounded  inland  by  highways  Fla.  92,  U.S.  Ul,  and  Fla.  29.  Flood 
Region  III  is  from  Everglades  City  to  Flaming  i,  Fla. , and  bounded 


39 


inland  by  highways  U.S.  Ul  and  Fla.  27.  This  region  is  the  major  part 
of  the  Everglades  National  Park.  The  last  flood  region,  IV,  is  the  re- 
maining portion  of  the  park  bounded  by  Florida  Bay  and  highways  Fla.  29 
and  U.F.  1. 

53.  Plate  39  compares  the  coastal  surge  envelope  for  a wall  and 
the  flooding  coast  boundary  conditions.  The  maximum  surge  for  both 
simulations  occurs  at  the  coast,  seaward  of  Everglades  City.  The  dif- 
ference between  the  maximum  surges  is  about  2 to  3 ft.  This  difference 
in  surge  envelopes,  however,  also  occurs  over  the  entire  length  of  the 
coast  boundary,  reflecting  a significant  flooding  condition.  In 
proper  perspective,  the  volume  of  water  flooding  the  coast  is  inconse- 
quential to  that  volume  entrained  in  the  nearshore  open  sea  circulation. 
Consequently,  water  can  be  removed  from  the  coast  and  be  almost  in- 
stantaneously replaced  without  significant  effect  on  the  nearshore  water 
level.  The  circulation  pattern,  however,  is  greatly  changed  due  to  the 
flooding.  The  water-surface  topography  together  with  the  depth-averaged 
velocity  field  for  a wall  and  the  flooding  coast  boundary  conditions 
are  presented  in  Plates  UO,  Ul,  and  h2.  In  the  early  surge  development 
(Plate  UO),  Florida  Bay  is  practically  drained  of  its  water.  As  the 
storm  center  passed  Cape  Sable,  the  surge  began  to  rapidly  develop  in 
Florida  Bay  (Plate  Ul).  The  effect  of  flooding  is  yet  to  be  noticed  in 
the  nearshore  water  levels  and  velocity  fields.  However,  Plate  U2 
which  is  a snapshot  6 hr  after  Plate  Ul,  clearly  shows  a very  different 
solution  in  the  nearshore  and  not-to-nearshore  conditions  resulting 
from  the  flooding  coast.  Of  note  is  the  enhancement  of  the  onshore 
velocity  component  which,  as  a result,  must  transport  more  sea  sediment 
to  (and  over)  the  coast  than  previously  computed  with  nonflooding  open- 
coast  models.  Furthermore,  as  the  mound  of  high  water  that  follows  the 
storm  upcoast  is  reduced  by  flooding,  that  effect  must  have  an  influence 
along  adjacent  nonflooding  coastal  segments.  This  occurs  at  Naples 
where  no  local  flooding  was  considered,  yet  flooding  that  occurred  many 
miles  away  results  in  a difference  of  0.5  ft  in  the  peak  surge  level. 
Certainly,  nonflooding  open-coast  models  are  Inappropriate  for  the 


study  of  Florida  Bay  or  the  Everglades  area.  One  should  also  question 
their  use  at  Naples. 

5^.  Table  9 compares  the  computed  water  level  envelope  and  that 
observed  for  various  coastal  locations.  The  water  level  envelope  is 
the  peak  water  level  from  the  time-record  made  by  the  superposition  of 
the  expected  tide  at  a coastal  location  and  the  surge  hydrograph  without 
regard  to  the  time  of  its  occurrence.  It  is  seen  that  the  model  results 
are  in  good  agreement  with  the  observed  water  level.  It  is  more  diffi- 
cult to  compare  computed  and  observed  inland  high-water  levels 

because  the  tide  cannot  be  simply  added  to  the  computed.  Furthermore, 
loses  some  meaning  in  surge  simulations  where  for  inland  locations, 
the  local  wind  setup  is  of  major  importance.  Consider  a surge  of  5 ft 
occurring  at  a time  of  low  tide  which  is  expected  to  be  -5  ft  relative 
to  msl.  The  surge  model  would  compute  flooding  relative  to  barrier 
elevations  based  on  msl  and  surge  levels  (rel  msl).  The  result  is  an 
erroneous  computation  for  and,  to  a lesser  degree,  also  for  the 

coastal  surge.  However,  for  small  tide  ranges  such  as  exist  in  the  Gulf 
of  Mexico,  the  computations  for  the  open-coast  water  le-v-l  are  repre- 
sentative of  the  superposition  of  tide  and  surge.  This  problem  can  be 
avoided  if  the  tide  function  is  input  at  the  sea  boundary  and  dynamically 
computed  along  with  the  surge.  It  is  recommended  that  this  should  be 
pursued  when  a generalized  predictive  tide  function  appropriate  to  the 
edge  of  the  continental  shelf  becomes  available.  Subject  to  the  above 
comments,  Table  10  shows  the  computed  (not  adjusted  for  the  tide)  and 
observed  high-water  levels  at  inland  locations.  The  observed  data  are 
taken  from  a poststorm  survey  which  are  presented  in  Plates  ^3--^5. 

Louisiana  Coast,  Hurricanes  Flossy,  Hilda,  Betsy, 
and  Carmen  Surge  Study 


55.  Hurricanes  Flossy,  Hilda,  Betsy,  and  Carmen  affected  the 
Louisiana  coast  from  Atchafalaya  Bay  to  the  Mississippi  River. 

Hurricane  Flossy  is  the  only  storm  in  this  study  that  paralleled  this 
coastal  region  with  subsequent  landfall  near  Pensacola,  Fla.  Tables  2-5 
present  the  storm  parameters  for  these  hurricanes.  Plate  k6  shows  the 


location  map,  the  four  storm  tracks,  and  the  computing  grid.  The 
Mississippi  River  delta  region  and  the  adjacent  levees  are  simulated  as 
nonovertopping  harriers  protruding  from  the  coast  in  the  manner  shown 
in  Plate  U6.  For  three  of  these  hurricanes,  Hydromet  has  estimated 
the  surface  wind  fields  at  various  storm  locations Information 
about  Hurricane  Carmen  is  from  personal  communication  witn  Dr.  Joe 
Pellissier,  National  Hurricane  Center  (NHC),  NOAA,  Miami,  Fla.  It  is 
noted  that  the  storm  track  as  observed  from  Hydromet,  other  weather 
bureau  offices  as  noted  in  poststorm  damage  survey  reports  of  the  U.  S. 
Army  Engineer  District,  New  Orleans  (NOD),  and  from  the  NHC  in-house 
historical  storm  track  record  are,  in  most  cases,  somewhat  different — 
in  particular,  the  track  of  Hurricane  Flossy  as  it  crossed  the 
Mississippi  River  delta  and  Hurricane  Carmen  as  it  approached  Lake 
Charles,  La.,  in  which  the  storm  was  decaying  (filling)  rapidly.  The 
storm  tracks  presented  in  this  report  are  from  the  NHC  record  except 
that  the  track  of  Hurricane  Hilda  prior  to  landfall  is  moved  6 n.mi. 
to  the  east  and  the  track  of  Hurricane  Flossy  in  the  Mississippi 
River  region  is  a compromise  between  the  NOD  report  and  NHC. 

56.  Six  discrete  ponding  areas  with  storage  area  information  and 
FHBC  data  are  presented  in  Table  11.  Flood  Region  I is  from  Marsh  Is. 
to  Point  Au  Fer  Is.  and  bounded  inland  by  elevated  railroads  or  highways, 
spoil  banks  of  canals,  bayous,  or  other  waterways.  Flood  Region  II 
is  from  the  terminal  point  of  Flood  Region  I to  Wine  Is.  and  bounded 
inland  by  the  spoil  banks  of  the  Hoimia  Navigation  Canal  to  the  west  and 
Bayou  Chene  to  the  east  and  highway  U.S.  90  to  the  north.  Flood 
Region  III  is  from  Wine  Is.  to  Belle  Pass  and  bounded  inland  by  the 
Houma  Navigation  Canal,  Bayou  Lafourche  to  the  east,  and  Bayou  Blue  and 
portion  of  U.S.  90  to  the  north.  Flood  Region  IV  is  from  Belle  Pass  to 
Barataria  Pass  and  bounded  inland  by  Bayou  Lafourche,  Barataria  Bay 
Waterway,  and  U.S.  90  to  the  north.  Flood  Region  V is  from  Barataria 
Pass  to  the  west  bank  of  the  Mississippi  River  levee  and  bounded  inland 
by  the  Barataria  Bay  Waterway  which  closes  with  the  levee  near  New 
Orleans,  La.  Flood  Region  VI  is  from  the  east  bank  of  the  Mississippi 
River  levee  to  the  entrance  of  Lake  Borgne  following  the  seaward  extent 
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of  the  marsh  area  and  bounded  inland  by  New  Orleans  and  U.S.  90.  All 
the  above-mentioned  bayous,  canals,  etc.,  have  spoil  banks  associated 
with  them  that  would  have  the  tendency  to  limit  flows  from  one  flood 
region  to  another. 

57.  For  Hurricane  Flossy  which  occurred  in  September  1956,  the 

coastal  surge  envelope  for  a wall  and  flooding  coast  boundary  condi- 
tions is  presented  in  Plate  UT.  The  comparison  of  computed  water 
level  envelope  and  observed  high-water  conditions  is  presented  in 
Table  12  for  coastal  areas.  Table  13  shows  the  peak  value  of  and 

the  observed  conditions  for  each  flood  region.  The  observed  conditions 
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are  taken  for  the  NOD  poststorm  survey  which  is  presented  in  Plate  i*8. 
The  comparison  of  computed  and  observed  conditions  is  not  as  good  as 
should  be  expected.  Although  the  surge  model  could  be  in  error,  more 
probably  it  is  either  the  winds  or  poor  field  survey  of  high-water 
conditions. 

58.  Hurricane  Hilda  occurred  in  October  I96U.  Plates  !+9-5^ 
show  the  water-surface  topography  together  with  the  depth-averaged 
velocity  field  for  a wall  and  flooding  coast  boundary  in  6-hr  intervals 
over  a 30-hr  period.  The  role  of  the  Mississippi  River  delta  in  limit- 
ing the  free  transport  in  the  continental  shelf  waters  is  clearly  seen 
in  the  series  of  the  above  snapshots.  The  effect  of  a flooding 
boundary  relative  to  a wall  is,  again,  seen  to  enhance  the  onshore- 
directed  transport  of  water.  The  principal  surge  area  occurs  at 
Caillou  Bay  as  seen  in  Plates  51  and  55,  which  show  the  surge  envelope. 
The  coastal  barrier  in  this  area  is  typically  a beach  berm,  generally 
less  than  3 ft  in  height  (msl).  Consequently,  significant  flooding 
occurred  from  I8OOZ  (Z  referring  to  Greenwich  Mean  Time),  3 October, 

to  approximately  0300Z  the  next  day.  The  winds  at  0300Z,  U October, 
are  directed  more  alongshore,  forcing  the  inland  surge  to  the  east. 

The  highest  observed  water  level  of  9.8  ft  occurred  inland,  to  the  east 
of  Caillou  Bay  along  the  Hoviraa  Navigation  Canal. 

59-  Tables  lU  and  15,  as  before,  compare  the  computed  and  ob- 
served high-water  levels  for  coastal  and  inland  regions,  respectively. 
The  observed  water  level  record  and  that  computed  (superposition  of 
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surge  and  expected  tide  as  determined  from  the  National  Ocean  Survey 

tide  table)  at  various  coastal  locations  are  presented  in  Plates  56-58. 

An  oscillation  of  8 hr  in  period  is  shown  in  the  computed  hydrograph  at 

Biloxi,  Miss.  (Plate  58).  The  period  of  oscillation  is  characteristic 

of  an  edge  wave.  Multioscillations  occur  from  the  trapping  of  energy 

between  the  right  lateral  boundary  and  the  protruding  Mississippi  River 

levee.  However,  other  computed  water  levels  in  this  region  do  not 

exhibit  such  oscillations  as  shown  in  Plates  57  and  59-  Plate  56  shows 

another  typical  problem  in  opeii-coast  surge  computations — that  of 

proper  wind  specification.  Before  storm  landfall,  the  analytically 

computed  winds  at  Eugene  Is.  are  always  directed  offshore,  resulting  in 

the  enhanced  drawdown  of  the  water.  Moreover,  the  island  is  on  the 

leeward  side  of  the  coast  and  land  frictional  influence  on  the  wind  is 

not  considered  in  the  hurricane  model.  Further  complication  occurs  in 

light  of  recent  developments  which  shows  that  the  surface  winds  in 

the  forward  section  of  the  storm  may  be  outflowing  (a  negative  ingress 
\ 2h 

angle).  The  results  presented  in  the  comparison  tables  and  in 
Plates  56-58  are  in  good  agreement  with  the  observed;  however,  improve- 
ment should  be  expected  with  a better  wind  model.  Plate  60  presents 
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the  observed  high-water  levels  as  reported  in  the  poststorm  survey. 

60.  Hurricane  Betsy  occurred  in  September  I965.  The  water- 
surface  topography  and  velocity  field  comparing  a wall  and  flooding 
coast  boundary  are  shown  in  Plates  6I-65.  The  effects  of  the  flooding 
boundary  on  the  nearshore  water  levels  and  currents  are  most  dramati- 
cally demonstrated  in  Plate  62.  Without  overtopping,  the  wall  condition 
is  trapping  water  in  excess  of  18  ft  at  the  corner  between  the  east 
Mississippi  levee  and  coastal  wall.  With  flooding  of  the  marsh  area, 
only  12  ft  of  water  is  computed  along  the  levee.  The  coastal  surge 
envelope  (Plate  66)  clearly  shows  the  effect  of  a flooding  boundary 
versus  that  with  a wall.  The  velocity  fields  are  considerably  differ- 
ent with,  again,  the  enhancement  of  the  onshore  velocity  for  the 
flooding  boundary.  Plate  62  also  shows  that  the  surge  is  Just  begin- 
ning to  develop  on  the  western  side  of  the  Mississippi  River  levee. 

Plate  63,  3 hr  later  than  that  snapshot  in  Plate  62,  shows  the  developed 


surge  in  this  area  while  to  the  east,  sustained  relatively  high  water  is 
continuing.  Plates  6h  and  65  show  the  decreasing  surge,  whereas 
f Plate  65  shows  the  reentrance  of  the  water  from  the  ponding  areas. 

61.  Tables  I6  and  IT  compare  the  computed  and  observed  high-water 
levels  for  coastal  and  inland  regions,  respectively.  Furthermore,  ob- 
served and  computed  hydrographs  at  two  coastal  locations  are  presented 
i 

1 in  Plates  67  and  68.  The  hydrograph  location  in  Plate  6?  is  to  the  left 

I of  landfall.  The  problems  associated  with  the  computed  large  drawdown 

have  been  discussed  previously.  No  presurge  water  level  anomaly  is 
I added  to  any  of  the  computed  resiilts  in  any  study  reported  herein, 

^ although  after  viewing  the  observed  hydrograph  at  Biloxi,  it  may  seem 

? appropriate.  The  results  from  this  surge  simulation,  more  than  in  any 

; other  presented  in  this  report,  show  the  necessity  in  treating  the  coast 

; as  a finite  barrier.  The  comparison  of  the  computed  and  observed  water 

i level  conditions  are  in  excellent  accord.  Plate  69  presents  the  observed 
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^ high-water  levels  as  reported  in  the  poststorm  survey. 

I 62.  Hurricane  Carmen  occurred  in  September  197^.  Snapshots  of 

the  water-surface  topography,  together  with  the  depth-averaged  velocity 
field  for  a wall  and  flooding  coast,  are  shown  in  Plates  70-72.  The 
principal  surge  developed  between  Caillou  Bay  and  the  entrance  to 
Terrebonne  Bay,  as  shown  in  the  surge  envelope  (Plate  73),  at  approxi- 
mately 09OOZ,  8 September  (between  the  time  of  Plates  70  and  81).  The 
^ surge  does  not  develop  in  Atchafalaya  Bay  until  I8OOZ.  After  this 

! time,  the  forward  motion  and  location  of  the  storm  are  in  question  due 

to  its  rapid  decay.  Since  the  surge  developed  in  East  Cote  Blanche 
Bay  during  the  decay  stage  of  the  storm,  it  is  viewed  as  fortuitous 
that  the  peak  computed  and  observed  water  level  are  in  good  agreement 
(Plate  7^).  The  phasing  of  the  peak  water  levels  at  Luke's  Landing, 

La. , is  in  poor  agreement  probably  due  to  the  storm  proceeding  faster 
on  the  track  than  that  used  in  the  simulation. 

63.  Plates  75-77  present  computed  and  observed  hydrographs  at 
other  coastal  locations.  The  observed  hydrographs  in  the  region  to 
the  east  of  the  Mississippi  River  (Plates  76  and  77)  show  a consider- 
able initial  water  elevation  at  the  start  of  computations.  This 
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elevation  is  not  totally  a presurge  condition  because  the  action  of  the 
wind  is  seen  to  reflect  the  observed  peak  levels  some  20  hr  later.  On 
the  other  hand,  this  initial  condition,  especially  at  the  Gulf  Outlet 
Canal,  is  not  due  to  the  hurricane  since  the  winds  on  the  continental 
shelf  are  very  weak  at  the  start  of  computations;  that  is,  if  the  surge 
simulation  had  begun  at  1200Z,  6 September  (2U  hr  in  advance  of  the 
starting  time  in  Plate  T6),  the  model  would  not  have  shown  over  h ft  of 
water  at  the  canal  entrance  at  1200Z,  7 September.  Other  processes 
are  occurring,  presumably  not  assignable  to  the  storm  surge,  for  what- 
ever caused  the  transient  water  level  bump  at  ITOOZ,  7 September 
(Plate  76)  is  probably  responsible  for  a similar,  well-correlated  bump 
at  2200Z,  7 September,  at  Biloxi  (Plate  77). 

64.  Tables  I8  and  19  compare  the  computed  and  observed  high- 
water  levels  for  coastal  and  inland  regions,  respectively.  There  are 
no  observations  along  the  coast  in  the  principal  surge  area  for  com- 
parison with  the  computed  value.  In  other  places  and  subject  to  the 
above  discussion,  the  comparison  between  computed  and  observed  water 
levels  is  in  good  agreement.  Plate  78  presents  the  observed  high-water 

27 

levels  as  reported  in  the  poststorm  survey. 
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PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 

65.  The  development  of  a two-dimensional,  time-dependent, 
curvilinear,  open-coast  storm  surge  and  tidal  model  is  presented.  In 
particular,  the  model  treats  the  coast  boundary  as  a finite  height 
barrier  which  is  broken  with  bay  entrances.  The  flooding  coast  routine 
in  its  treatment  of  submerged  and  exposed  (on  one  side)  barriers  is 
considered  conservative.  The  results  from  the  model  compare  most 
favorably  with  those  obtained  from  two  hydraulic  physical  models  in- 
volving tidal  or  pseudosurge  simulations.  These  experiments  revealed 
that  the  model  could  be  improved  in  its  treatment  of  open  sea,  finite 
height  structures  such  as  weirs  or  jetties. 

66.  The  results  of  five  storm  surge  simulations  are  presented  in 
which  observed  and  computed  conditions  are  compared.  These  simulations 
employed  two  separate  stretches  of  coastline,  each  of  different  coastal 
characteristics.  Again,  the  model  results  compare  most  favorably  with 
those  ooserved.  However,  from  currently  available  data,  a relatively 
low  value  for  the  wind  drag  coefficient  at  high  wind  speeds  is  used. 

This  is  necessary  to  compensate  for  the  liberal  estimate  for  the  maximum 
winds.  In  hindcasting  situations,  it  is  not  guaranteed  that  the  same 
procedure  is  used  for  estimating  the  maximimi  wind.  The  wind  model 
employed  in  these  surge  simulations  from  historical  storms  requires 
input  of  the  maximum  wind.  The  problem  can  be  resolved  by  using  a wind 
model  based  on  the  physics  of  the  marine  boundary  layer  that  will 
provide  a wind  field  as  a result  of  dynamic  meteorological  computations. 
It  is  recommended  that  such  a model  be  obtained  (as  they  do  exist)  and 
incorporated  into  the  hydrodynamic  flooding  coast  model.  The  resulting 
composite  model,  therefore,  would  be  appropriate  for  both  lorecasting 
and  hindcasting  and  involve  a drag  coefficient  reflective  of  the  cur- 
rently available  data. 

67.  The  following  recommendations  are  secondary  to  that  of  in- 
corporating a dynamic  wind  model  in  SSURGE  III.  It  can  be  shown  that 
the  astronomical  tide  does  distort  the  predicted  coastal  and  inland 
flood  water  levels  when  it  is  not  incorporated  into  the  dynamic 


computations;  that  is,  equivalent  results  are  not  obtained  by  the  super- 
position of  surge  and  tide  and  that  simulation  where  the  tide  is  in- 
cluded in  the  sea  boundary  forcing  condition.  For  the  east  coast  of 
the  United  States  where  there  are  large  tidal  ranges,  the  tide  should 
be  included  in  the  computation.  The  results  of  the  five  storm  surge 
simulations  are  considered  sufficiently  accurate  using  superposition  as 
they  occurred  in  Gulf  of  Mexico  waters  where  the  tidal  range  is  small. 

It  is  recommended  that  a generalized  tide  prediction  function  which  is 


appropriate  to  the  edge  of  the  continental  shelf  be  incorporated  into 
the  surge  computations. 

68.  It  is  recognized  that  other  approaches,  equations,  and  co- 
efficients are  reported  or  can  be  developed  for  the  treatment  of  a 
flooding,  finite  height  barrier  coast.  These  should  be  investigated 
and  the  more  promising  approach  tested  against  the  one  reported  herein 
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Table  7 

^LTICE  in  Finite  Hcluht  CoMt  Boundarv_D«u  for  Murrell,  Inlet 


Table  9 


Comparison  of  Computed  and  Observed  Coastal 
High-Water  Levels  for  Hurricane  Donna 


Computed* 

Observed 

Coastal 

High-Water 

High-Water 

Location 

Level  (ft) 

Level  (ft) 

Average 

Estero  Bay 

8.5 

7.4,  7.7,  7.8, 

8.1,  9.5 

8.1 

10.5 

9.7,  10.4 

10.1 

11.0 

9.8,  10.2, 

10.4,  11.1 

10.4 

8.3  n. mi.  North 
of  Naples,  FL 

11.2 

8.1,  10.2, 

11.9 

10.1 

3.3  n.  mi.  North 
of  Naples,  FL 

10.2 

10.0 

10.0 

Naples,  FL 

9.8 

(w/o  flood- 
ing coast, 
10.4) 

8.1,  9.4,  10.0, 
11.3,  11.6 

10.1 

Marco  I. 

(1.2  n.mi.  inland) 

9.5 

10.2 

10.2 

3.8  n.mi.  East 
of  Cape  Romano 
(1.5  n.mi.  inland) 

11.7 

7.8,  11.7 

9.8 

Everglades  City,  FL 

15.0 

N/A 

Shark  River 

9.3 

N/A 

Flamingo,  FL,  on 

Cape  Sable 

11.0 

(w/o  flood- 
ing coast, 
12.5) 

12.0 

12.0 

(Continued) 


Adjusted  for  astronomical  tide 


Table  9 (Concluded) 


Coastal 

Location 

Computed* 
High-Water 
Level  (ft) 

Observed 

High-Water 

Level  (ft) 

Average 

Florida  Keys 

Rock  Harbor,  FL 

10.1 

5.6, 

7.0,  8.9 

7.2 

Plantation  Key 

Upper  Mat.  Key 

8.2 

4.0, 

9.0, 

4.8,  6.3, 
10.0,  10.6 

7.5 

Craig,  FL 

7.3 

8.5 

8.5 

Long  Key 

6.4 

6.0 

6.0 

Grassy  Key 

8.3 

7.0, 

7.4 

7.2 

Vaca  Key 

9.9 

7.5, 

8.5, 

o 

00 

O 00 

00  00 

8.2 

Torch  Key 

7.0 

3.4 

3.4 

Sugarloaf  Key 

5.1 

6.6 

6.6 

10.0  n.mi.  East 
of  Key  West,  FL 

4.4 

4.9 

4.9 

5.5  n.mi.  East 
of  Key  West,  FL 

3.5 

4.7 

4.7 

Key  West,  FL 

2.3 

1.9 

1.9 

L 


* Adjusted  for  astronomical  tide 


Table  10 

Comparison  of  Computed  and  Observed  Inland 
High-Water  Levels  for  Hurricane  Donna 


Inland 

Location 

Computed* 
Ponding  Region 
Level  (ft) 

Observed 
High-Water 
Level  (ft) 

Marco  Is. 

(1.2  n.mi.  inland). 

Flood  Region  I 

5.9 

5.3 

(7.7  n.mi.  inland 
near  FL  41) 

Everglades  City,  FL 

Flood  Region  II 

7.1 

6.2 

(7.0  n.mi. 
8.5,  9.7 
(5.1  n.mi. 

inland) 

inland) 

Shark  River, 

Flood  Region  III 

6. 4 

6-6.5 
(4.2  n.mi. 
and  0.5  n. 
river) 

inland 
mi.  from 

Flamingo,  FL,  on 

Cape  Sable, 

Flood  Region  IV 

7.0 

N/A 

j 


I 

I 


* Not  adjusted  for  astronomical  tide. 


Table  12 

Comparison  of  Computed  and  Observed  Coastal 
High-Water  Levels  for  Hurricane  Flossy 


Coastal 

Location 

Computed* 

High-Water 

Level  (ft) 

Observed 
High-Water 
Level  (ft) 

6.1  n.mi.  East  of 

East  Timbalier  1. 
at  Bay  Champagne 

2.3 

3.7  (EST**) 

Southwest  Pass  Ent. 

2.0 

3.7 

Quarantine  Bay 

9.6 

10.8,  11.9,  12.1 

Pointe  A La  Hache,  LA  to 

Point  Pleasant,  LA 

8-9 

10.3,  10.5  (EST) 
'10.7,  11.0  (EST) 

Breton  I. 

6.0 

7 to  8 (EST) 

* Adjusted  for  astronomical  tide. 

**  Furnished  by  various  oil  companies  and  by  the  Freeport  Sulphur  Co. 


Table  13 

Comparison  of  Computed  and  Observed  Inland 
High-Water  Levels  for  Hurricane  Flossy 


Computed* 

Inland  High-Water 

Location  Level  (ft) 


Observed 
High-Water 
Level  (ft) 


Barataria  Bay  to  0.4  3.6**,  3.6**,  5.2**,  5.5**, 

Mississippi  River  Levee,  6.9**,  7.3**,  8.0** 

Flood  Region  5 


Lake  Borgne, 
Flood  Region  6 


2.5  4.0,  5.2,  6.3,  8.5,  6.2t,  6.4t, 

6.9t,  7.7t,  8. It,  8.6t,  9.4t, 
10. 9t 


; 

■ 

i 

i 


I 

i 


1 

\ 

•I 

I 


NOTE:  Flood  Regions  1-4  not  affected.  i 

* Not  adjusted  for  astronomical  tide.  j 

**  Observations  were  located  along  landward  side  of  coastal  j 

barrier.  j 

t Observations  were  located  along  southwest  side  of  Lake  Borgne.  j 


r 


Table  14 

Comparison  of  Computed  and  Observed  Coastal 
High-Water  Levels  for  Hurricane  Hilda 


Coastal 

Location 

Computed* 

High-Water 

Level  (ft) 

Observed 
High-Water 
Level  (ft) 

Gulf  Ent. , Freshwater  Bayou 

1.5 

2.3 

East  Cote  Blanche  Bay  at 

2.5 

4.4 

Lukes  Landing,  LA 

Atchafalaya  Bay  at  Eugene  I. 

3.8 

3.3 

(w/o 

flooding  coast,  4. 

3) 

Lake  Pelto,  2.4  n.mi.  Landward 

6.8 

7.4 

from  Ent.  to  Terrebonne  Bay 

(w/o 

flooding  coast,  7. 

2) 

Grand  I. 

3.2 

4.0 

(w/o 

flooding  coast,  3. 

2) 

East  Side  Barataria  Bay  near 

5.0 

5.5 

Gulf  Ent. 

Southwest  Pass,  Mississippi 

1.5 

3.2,  3.5 

River 

Quarantine  Bay  at  Ostrica,  LA 

4.0 

4.6 

(w/o 

flooding  coast,  4, 

4) 

Ent.  Mississippi  River  - Gulf 

4.0 

4.6 

Outlet  Canal 

(w/o 

flooding  coast,  4, 

4) 

Biloxi,  MS 

5.2 

4. 6 

(w/o 

flooding  coast,  5. 

.2) 

* Adjusted  for  astronomical  tide. 


k i 


Table  15 

Comparison  of  Computed  and  Observed  Inland 


High-Water  Levels  for  Hurricane  Hilda 


Computed* 

Observed 

Inland 

High-Water 

High-Water 

Location 

Level  (ft) 

Level  (ft) 

Atchafalaya  Bay, 

Flood  Region  I 

0. 

3.8, 

5.1 

Atchafalaya  River 
to  Houma  Canal, 

Flood  Region  II 

3.8 

None 

Terrebonne-Timbalier  Bay, 

3.4 

3.0, 

3.8,  4.0, 

4.7, 

Flood  Region  III 

5.4, 

6.6,  6.9, 

6.9, 

7.0, 

7.4,  9.8 

Bayou  Lafourche 
to  Baratarla  Bay, 

Flood  Region  IV 

0.5 

None 

Barafaria  Bay  to 

Mississippi  River  Levee, 

Flood  Region  V 

1.6 

4.0 

Lake  Borgne, 

3.3 

3.8, 

4.5,  4.7, 

5.0, 

Flood  Region  VI 

5.2, 

5.3 

* Not  adjusted  for  astronomical  tide. 


Table  16 


Comparison  of  Computed 


and  Observed  Coastal 


High-Water  Levels  for  Hurricane  Bets:, 


Coastal 

Location 


Computed* 
High-Water 
Level  (ft) 


East  Cote  Blanche  Bay  at 
Lukes  Landing,  LA 


4.5 


Observed 
High-Water 
Level  (ft) 


3,8 


Grand  I. 


7.2 

(w/o  flooding  7.4) 


8.8 


East  Side  Barataria  Bay  o-J 

near  Gulf  Ent.  (w/o  flooding  coast,  9.8) 


West  Side  Mississippi 
River  Levee;  Empire,  LA, 
to  Venice,  LA 


6.2  to  12.2  5.7,  7.4,  7.7, 

(w/o  flooding  coast,  8.8,  9.2,  10.4 
6.2  to  16.1) 


East  Side  Mississippi 
River  Levee;  Pointe  A La 
Hache,  LA,  to  Ostrica 
(Brenton  Sound) , LA 


12.4  to  15.2  13.6,  13.7, 

(w/o  flooding  coast,  14.4,  14.5,  15.7 

12.8  to  23.1) 


Gulfport,  MS 


9.0  (Approx.)  10.7 


Biloxi,  MS 


7.6 


8.6 


* Adjusted  for  astronomical  tide. 


Table  17 


Comparison  of  Computed  and  Observed  Inland 
High-Water  Levels  for  Hurricane  Betsy 


Inland 

Location 

Computed* 
High-Water 
Level  (ft) 

Observed 

High-Water 

Level  (ft) 

Atchafalaya  Bay, 

Flood  Region  I 

0. 

2.9,  3.9 

Atchafalaya  River 
to  Houma  Canal , 

Flood  Region  11 

0.9 

3.0,  4.3 

(along  Houma  Canal) 

Terrebonne-Timbalier  Bay, 
Flood  Region  III 

1.8 

2.8,  3.0,  4.3 

Bayou  Lafourche 
to  Barataria  Bay, 

Flood  Region  IV 

0.5 

3.4,  5.4 

Barataria  Bay  to 

Mississippi  River  Levee, 
Flood  Region  V 

2.1 

3.4,  5.0,  5.3 

Lake  Borgne, 

Flood  Region  VI 

7.1 

6.4,  8.3,  8.8,  9.1, 
9.3,  9.8,  10.1,  10.7, 
11.0,  11.7,  14.4 

* Not  adjusted  for  astronomical  tide. 
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Table  18 

Comparison  of  Computed  and  Observed  Coastal 


Coastal 

Location 

Com.puted* 
High-Water 
Level  (ft) 

Observed 

High-Water 

Level  (ft) 

East  Cote  Blanche  Bay  at 
Lukes  Landing,  LA 

4.3** 

4.5 

Atchafalaya  Bay  at 

Lower  Atchafalaya  River 

4.5** 

4.9 

Atchafalaya  Bay  at 

Eugene  I. 

3.8 

(w/o  flooding  coast. 

4.0) 

3.7 

Grand  I. 

4.7 

(w/o  flooding  coast. 

4.8) 

4.2 

Ent.  Mississippi  River  - 
Gulf  Outlet  Canal 

5.1 

(w/o  flooding  coast. 

5.4) 

5.7 

Biloxi,  MS  4.0 

(w/o  flooding  coast,  no 

change) 

4.5 

Caillou  Bay 

9.2 

(w/o  flooding  coast. 

11.6) 

None 

Ent.  Terrebonne  Bay 

8.7 

(w/o  flooding  coast. 

11.1) 

None 

1 

1 

] 

j 

* Adjusted  for  astronomical  tide. 

**  The  time  of  high  water  occurred  after  1200  CST  8 Sept  1974.  The 
water  level  is  approximate  because  the  structure  of  the  storm  after 
1200  CST  is  uncertain  due  to  its  rapid  decay.  Computations  ended 
at  1800  CST  8 Sept. 


Comparison  of  Computed  and  Observed  Inland 
High-Water  Levels  for  Hurricane  Carmen 


Inland 

Location 

Computed* 

High-Water 

Level  (ft) 

Observed 
High-Water 
Level  (ft) 

Atchafalaya  Bay, 

Flood  Region  I 

0.6** 

2.8,  3.2,  3.8, 
4.0,  4.8 

Atchafalaya  River  to 

Houma  Canal, 

Flood  Region  II 

3.8 

3.3 

Terrebonne-Timbalier  Bay, 

Flood  Region  III 

4.4 

2.9,  3.8,  4.6, 

5.9,  11.6 

Bayou  Lafourche  to 

Barataria  Bay, 

Flood  Region  IV 

0.5 

None 

Barataria  Bay  to 

Mississippi  River  Levee, 

Flood  Region  V 

1.4 

3.1 

Lake  Borgne, 

3.4 

5.0,  6.0 

Flood  Region  VI 


* Not  adjusted  for  astronomical  tide. 

**  Computed  bay  level  rising  at  end  of  computation.  See  footnote  ** , 
Table  18. 


COMPUTING  GRID  AND 
LOCATION  MAP 


GALVESTON  entrance  CHANNEL  AREA  WITH 
NUMERICAL  PROTRAYAL  OF 
ENTRANCE  CHANNEL  JETTIES 
time  : 4-4  hours 


velocity,  fps 


30  40 

Time,  hr 


GRID  POINT  (M.  «),  GALVESTON 

entrance  channel 

GRID  POINT  (16,  9X  rollover  PASS 


NUMERICAL  MODEL  DEPTH 

averaged  onshore 
velocity  component 


GALVESTON  BAY  MODFL 

HURRICANE  CARLA 
SIMUL  ATION 
VELOCITY 


PLATE  9 


Time  , HR 


LEGEND 

GRID  POINT  (II,  9), GALVESTON 

entrance  channel 

GRiO  POINT  (16,9),  ROLLOVER  PASS 


numerical  model  depth 
averaged  onshore 
velocity  component 


GALVESTON  bay  model  STUDY 


LARGE  RADIUS, 

HIGH  TRANSLATION 
HURRICANE  SIMULATION 
VELOCITY 


HID-AUH/  AKMT  tNblNEER  WATERWATS  EXPERIMENT  STATION  VICKSBUR6  MISS  F/G  4/2 

^ NEARSHORE  NUMERICAL  STORM  SURGE  AND  TIDAL  SIMULATION. (U) 

SEP  77  J J WANSTRATH 

unclassified  WES-TR-H-77-17  NL 


MURRELLS  INLET  MODEL  STUDY 

9-FT  SURGE  SIMULATION 
HYOROGRAPH  COMPARISON 
MURRELLS  INLET  ENTRANCE 


TIME,  HR 


LEGEND 

MID  depth  velocity  \ physical 

AT  STATION  1 / MODEL 

DEPTH  AVERAGED  ON  Y , ^ 

SHORE  VELOCITY  AT  ( NUMERICAL 
GRID  POINT  (7,  7)  J MODEL 


MURRELLS  INLET  MODEL  STUDY 

9-FT  SURGE  SIMULATION 

VELOCITY  comparison 


14 


LEGEND 

— H AT  GRID  POINT  (7,  1) 
^10 
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H = Hr 
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Hr  = WATER  ELEVATION  DUE  TO  RADIATION  COMPONENT 


time  5 10  HR: 


HlO  = ^0  E 

Ksl 


A=  AMPLITUDE,  FT 
T=  PERIOD,  HR 

t = time 

0=  PHASE,  DECREES 


Ak 

l»K 

5.26 

- 

- 

6.12 

10.0 

-177 

0.98 

5.0 

-20 

0.30 

3.3 

114 

0.05 

2.5 

35 

time  > 10  HR,  Hio=0 
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NUMERICAL  MODEL  SEA 
BOUNDARY  FORCING  FUNCTION 


PLATE  26 


LEGEND 


HYOROGRAPH 

OPEN  COAST 


VELOCITY.  FPS 


LEGEND 

AT  STATrON''^''°‘^''''']'’”"SICAL  MODEL 

- — DEPTH-AVERAGED  I 

ON-SHORE  VELOCITY  InomERICAL  MODEL 
AT  GRID  POINT  (7,  7)J 


MURRELS  INLET  MODEL  STUDY 

12-FT  SURGE  SIMULATION 

VELOCITY  COMPARISON 
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APPENDIX  A:  NOTATION 
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Storage  area  of  the  attached  bay  known  as  a function, 

Channel  cross-sectional  area  at  mean  sea  level 

Surface  area  at  mean  sea  level  of  the  grid  cell  representative 
of  H 

Surface  area  of  the  attached  bay  at  mean  sea  level 
Transformation  coefficients 

Nondimensional  variable  wind  drag  coefficient 
Transformation  coefficient 

Nondimensional  channel  discharge  coefficient 

Nondimensional  coefficient  for  an  overtopping  barrier  which  is 
exposed  on  one  side 

Nondimensional  coefficient  for  a submerged  barrier 
Depth  of  water  (H  - D^) 

Local  water  depth  relative  to  mean  sea  level 
Coriolis  parameter 

Scale  factor  associated  with  the  orthogonal  curvilinear  co- 
ordinate system 

Acceleration  due  to  gravity 

Water  level  relative  to  mean  sea  level 

Water  level  relative  to  mean  sea  level  in  the  ponding  area 

Hydrostatic  elevation  of  the  sea  surface  corresponding  to  the 
departure  of  the  atmospheric  pressure  from  ambient 

Predicted  water  level  at  the  coast  barrier  without  correction 
for  flooding 

Number  of  grid  points  per  computational  lattice  along  the  S 
and  T axes 

Nondimensional  variable  seabed  drag  coefficient  that  depends 
on  the  seabed  condition  and  water  depth 

Length  of  coastline  at  unit  elevations,  ^ ~ 0,1,2... 

Map  factor  relating  prototype  length  to  x,y  units 
Atmospheric  pressure 

Central  atmospheric  pressure  of  hurricane 
Far-field  atmospheric  pressure 


A1 
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S ,S 
n’  p 


T ,T 
x’  y 


W ,W 
x’  y 


W, 


R 


W. 


QgiQip  Voliome  transport  per  unit  width  in  the  S and  T directions  or, 
equivalently,  in  the  local  direction  of  C and  n in  proto- 
type space 

r„  Distance  from  hurricane  center  to  H grid  points 
rl 

Distance  from  the  storm  center  to  the  region  of  maximum  winds 
Ordinate  axis  of  the  stretched  shelf  coordinate  system 
Distance  normal  to  the  seaward  hourdary  and  along  the  coast 
Abscissa  axis  of  the  stretch  shelf  coordinate  system 
Long-wave  travel  time 

The  X and  y components  of  the  forward  speed  of  this  hurricane 
center 

Bay  entrance  channel  width 

The  X and  y .wind-speed  components  for  a moving  hurricane 
The  maximum  wind  speed  for  the  stationary  hurricane 
Wind  speed  at  an  elevation  10  meters  above  the  water  surface 
Rectilinear  coordinate  in  prototype  space 

The  X coordinate  of  the  hurricane  center  in  prototype  space 
The  X coordinate  of  the  H grid  points  in  prototype  space 
Rectilinear  coordinate  in  prototype  space 

The  y coordinate  of  the  hurricane  center  in  prototype  space 
Yp.  The  y coordinate  of  the  H grid  points  in  prototype  space 

O 

Z,  Unit  elevations  of  coastal  barrier  as  denoted  by  k = 0,1,2... 
b 

6 The  extent  of  q ( + ) in  the  image  space  of  (C,ti) 

AS  Constant  spaced  computing  grid  increment  between  dependent 
variables  in  the  S direction 

At  Aigorithm  time  increment  between  computing  lattices 

AT  Constant  spaced  computing  grid  increment  between  dependent 
variable  in  the  T direction 

n The  offshore  curvilinear  coordinate  in  prototype  space  and  the 
abscissa  in  the  image  space  of  (C,n) 

0 The  local  angle  between  the  C and  x axis  in  prototype  space 

X The  positive  V.ori zontal  extent  of  the  region  to  be  mapped  in 

prototype  space  and  constrained  to  equal  the  extent  of  f; 
in  the  image  space. 

U,v  Functions  transforming  C and  n to  the  stretched  shelf 
coordinate  system  (S,T) 
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